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O 

' For uniform systems of spin-less fermions in d spatial dimensions with d > 1, interacting through 

^SJ ' the isotropic two-body potential v{r — r'), a celebrated theorem due to Luttinger (1961) states that 

> I , under the assumption of validity of the many-body perturbation theory the self-energy S(fc;e), with 

C<1h' < k^Skp (where kp stands for the Fermi wavenumber), satisfies the following universal asymptotic 

, relation as e approaches the Fermi energy ep'- ImE(fc; e) ~ =pafc(e — ef)'^, e^sf, with Uk > 0. As this 

is, by definition, specific to self-energies of Landau Fermi-liquid systems, treatment of non-Fermi-liquid 
systems are therefore thought to lie outside the domain of applicability of the many-body perturbation 
theory; that, for these systems, the many-body perturbation theory should necessarily break down. We 
demonstrate that ImE(fc; e) ~ =pQfe(e — ef)^, e^ef, is implicit in Luttinger's proof and that, for d > 1, 
^ , in principle nothing prohibits a non-Fermi-liquid-type (and, in particular Luttinger-liquid-type) E(fc;e) 

from being obtained within the framework of the many-body perturbation theory. We in addition indicate 
how seemingly innocuous Taylor expansions of the self-energy with respect to k, e or both amount to 
^ , tacitly assuming that the metallic system under consideration is a Fermi liquid, whether the self-energy 

is calculated perturbatively or otherwise. Proofs that a certain metallic system is in a Fermi-liquid state, 
based on such expansions, are therefore tautologies. 
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I. INTRODUCTION 



> 

\l ' Landau's (1957) phenomenological theory for the low-lying excited states of many-particle systems is based on the 
assumption that these states are characterised by a distribution of a dilute gas of quasi-particles (QPs), implying 
that interaction amongst these QPs may, to a good approximation, be neglected (Nozieres 1964, Pines and Nozieres 
1966, Abrikosov, Gorkov and Dzyaloshinski 1975) (for a comprehensive review see Platzman and Wolff (1973)). [] 
This ideal situation corresponds to the case in which QPs are the exact one-particle eigenstates of the many-body 
^ Hamiltonian, or, equivalently, the energies of the QPs are infinitely sharply defined. The success of Sommcrfeld's 
independent-particle model (for example Ashcroft and Mermin (1981, Ch. 2)) for conduction electrons in simple 
metals, despite the apparent non-negligible strength of the electron-electron Coulomb interaction, can be viewed as 
a most strong pillar of Landau's phenomenological theory. Incorporation of the residual two-body interaction among 
Q , the ideal QPs (at least beyond the Hartree-Fock or the static exchange approximation), which evidently should be 
^ ' weaker than the bare Coulomb interaction (insofar as low-energy scattering processes are concerned), ^ renders the 
• ^ , QP states non-stationary. In other words, the actual QPs, if such entities can at all be meaningfully defined (§ III), 
' are merely approximate one-particle eigenstates of the many-body Hamiltonian; the non-stationary nature of a QP 
^ , 'eigenstate' signifies that its energy, in contrast with that of an exact eigenstate, is un-sharp, or the spectral function 
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^ To be precise. Landau's Fermi-liquid theory is applicable for temperatures T in the range {Tc, T*), where T* and Tc stand, 
respectively, for the coherence and the transition temperature; for T below Tc the system is in a superconducting state (Kohn 
and Luttinger 1965). For details as well as a comprehensive discussion of instabilities of Fermi liquids see (Metzner, Castellani 
and Di Castro 1998). 

^ Otherwise the independent-particle model would not be an appropriate zeroth-order approximation. 
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of QPs consists of peaks with finite widths. Thus the true, that is interacting, QPs correspond to superpositions of 
nearly degenerate (as well as degenerate) many-body eigenstates, so that an initially well-defined QP ceases to be 
particle-like (in the sense of possessing a reasonably well-defined energy) with the passage of time. 

The information with regard to QPs is contained in the single-particle Green function G. The formal Lehmann 
(1954) spectral representation of this function in the energy domain reveals that energies of the single-particle ex- 
citations, the QPs, can be identified with the singular points (often inappropriately designated as 'poles' — § V) 
along the real energy axis (e-axis) (Fetter and Walecka 1971). Such observation is of little practical relevance when 
systems in the thermodynamic limit are concerned. For macroscopic systems, the continuum of the energy levels 
^ gives rise to branch points (which are not isolated singularities) and branch cuts along the energy axis. This is 
naturally consistent with the above-indicated observation with regard to the broad rather than delta-function-like 
spectral function of QPs. Such a dramatic change in the analytic structure of G(e) gives rise to a number of effects 
the subtleties of which are often overlooked in the literature. In this work we discuss in some detail a number of these, 
insofar as they are relevant to the main objective of our work. One of these effects is associated with the fact that 
whereas a reasonably sharp peak in the spectral function may be described by a single pole with some small imaginary 
part, such description is in violation of the fact that G(e) is a perfectly regular function on the entire complex energy 
plane, excluding the real e-axis. In § III we show that such complex poles do not simply correspond to the analytic 
continuation of G{s) into the complex energy plane (i.e. the physical Riemann sheet (RS)), rather to the analytic 
continuation of the latter to a RS; we refer to this as a non-physical RS (see Appendix A). 

The main complication, from the restricted viewpoint of the present paper, arises from the branch points of the 
Green function G{k] e), in particular that ed e — ep which can be shown to be also a branch point of the self-energy, 
S(fc;e). From this it follows that, for instance, such expression as (Hugenholtz 1957, DuBois 1959b, Luttinger 1961) 
ImS](fc;£) ^ =pafc(£ — Sp)"^, £^£_f, cannot be considered as the leading term in the Taylor series expansion of the 
self-energy operator at e — £p. Thus there exists no a priori reason for the dispersion of the QP energies to be a 
'smooth' function of fc in a neighbourhood]^ of kp (for specification see further on) similar to that of non- interacting 
QPs, namely 0£j! oc fc^, with fc:=||k||. In fact, for the special case of 'Luttinger liquids' (Haldane 1980, 1981, Luttinger 
1963, Mattis and Lieb 1965, Dover 1968, Anderson 1997), without hereby specifying the dimension of the spatial space 
to which the corresponding system is confined, or marginal Fermi liquids (Varma, et al, 1989, Littlewood and Varma 
1991, Kotliar, et al., 1991), we explicitly demonstrate that £k is not differentiable at k = kp (see §§ IV. C and VI). As 
we shall emphasise later in this paper (§ IV. C), QPs in Fermi liquids are special in that their energy dispersion in the 
close vicinity of the Fermi surface is, apart from a scaling factor arising from the re-normalisation of the electron mass 
or the Fermi velocity, non-interacting- like. Our close inspection of the proof of the Luttinger (1961) theorem (§ IV), the 
contents of which we have spelled out in the Abstract, reveals that this 'non-interacting-like' assumption with regard to 
the dispersion of the QP energies is implicit in Luttinger's (1961) proof, so that Luttinger's theorem in essence amounts 
to a statement concerning the consistency of this assumption with the property ImS(fc;e) ~ =pafc(£ — Ef)^, e^ef, 
for s —>■ Sp. In view of our rigorous demonstration with regard to non-differentiability oi Sk at k = kp for Luttinger 
liquids, or marginal Fermi liquids, which invalidates Luttinger's implicit assumption, we are left to conclude that in 
principle nothing precludes existence of a non-Fermi-liquid behaviour for metallic systems f\ in spatial dimensions 



^ From dimensional considerations it follows that for a system with linear dimension L, the separations between the energy 
levels scale like (Landau and Lifshitz 1980, p. 14). 

* At places in the present text, such as here, we refer to 'QP energies' even though QPs may not be well-defined. This is 
justified by the fact that even in these cases the equation for the 'QP energies' (i.e. Eq. (IL4| below) does have a solution — the 
ill-defined nature of QPs in these cases is associated with the fact that the single-particle eigenstates corresponding to these 
energies (see Eq. ( |lo[ ) below) do not describe particle-like excitations (for details see in particular § IV. C). 

^ Throughout the present work by 'neighbourhood' we refer to an open non-vanishing interval (e.g. along the real axis) or 
domain (on the complex plane) which, however, may be arbitrarily small. 

® We have e5!:=ft^fc^/[2me] = e'l^ + [h^kp /mg){k — kp) + {Ti? /\lm^){k — kp)'^. The important feature of this energy dispersion 
is that it is of the general form el = + <;{k ~ kp) + oik — kp) for k ^ kp (see § IV. C), where <; is a finite constant (explicitly, 
in the case at hand we have ? = h^kp/nie) and o{k — kp) stands for a function that in comparison with (fc — kp) is vanishingly 
small as k ^ kp. In this respect the property oc k^ is to us of relatively minor relevance. We point out that the coefficient 
of the linear term in the QP energy dispersion (here ?), if such a coefficient at all exists, has the relevance that through it such 
attributes as the Fermi velocity vp and effective mass m* can be assigned to the corresponding QP (see § IV. C, specifically 
Eqs. (|^ and ^). 

Our emphasise throughout this work on the metallic nature of the systems under consideration corresponds to the fact that 
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larger than unity. In particular, contrary to the general belief, the break-down of the many-body perturbation theory 
(Luttinger 1961, Anderson 1988, 1989, Varma, et al. 1989, Anderson 1990a,b, 1991, Littlewood and Varma 1991, 
Kotliar, et al. 1991, Anderson 1992, 1993, 1997) for a metallic system and the nature of the low-lying single-particle 
excitations in this stand in no direct relationship. ^ 

For clarity of discussions, we briefly specify the scheme by which we classify a metallic system as a Fermi liquid. 
This scheme, which to our best knowledge has not been utilised earlier in a similar context, has the advantage of 
being general and independent of the details of the metallic system to which it is applied. To this end, we first point 
out that the behaviour Iml](/c;e) ~ To^kis — £f)^, e^ef, for e ep, which is specific to Fermi liquids, leads to (see 
§§ IV. C, V and Appendix C) ReS(fc,e) — S(A;,ei?) ^ /3fe(e — ep) for e ep, with (3k a constant to be specified in 
Appendix C. The distinguishing aspect of I](fc,e) whose real and imaginary parts have such asymptotic forms is as 
follows. 

Condition A: E(/ci?;e) is a continuously differentiable function ^ of e in a neighbourhood of ep. 

This condition suffices to guarantee a finite discontinuity Zkp in the momentum-distribution function n(fc) at k = kp 
(§ V). A finite Zkp is generally considered to signify a Fermi liquid. If this merely evidences a Fermi liquid, then 
condition (A) would constitute the necessary and sufficient condition for a metallic system to be one (see § V). However, 
condition (A) can be shown (see §§ IV. C and V) not to be sufficient for the QP energy to possess a 'well-defined' 
dispersion in a neighbourhood of A; — kp, that is one which is expressible in terms of a continuously differentiable 
function of k in this neighbourhood. For the existence of such a dispersion, the following condition is also required to 
be satisfied. 

Condition B: Y,{k;ep) is a continuously differentiable function of A; in a neighbourhood of fc = kp. 

We therefore choose to consider a metallic system as a Fermi liquid provided the corresponding E(fc;e) satisfies 
both condition (A) and condition (B); this specifies the notion of 'smooth' introduced above. According to this, which 
is physically the most sound definition for Fermi liquids, ^ a finite Zkp, although necessary, is not sufEcient for a 
system to be a Fermi liquid. We explicitly demonstrate (§ VI) that whereas T,{k;ep) corresponding to a Fermi liquid 
is (by definition) a continuously differentiable function of fc in a neighbourhood of k = kp, the same may not apply 
to I](fc;e) when e ^ ep. We further demonstrate (§ VI) that S(fc;e) pertaining to a metallic system (whether Fermi 
liquid or otherwise), which is not continuously differentiable with respect to fc in a neighbourhood of fc = fco (where 
fco would possibly be fc^) for some e ^ ep, is not continuously differentiable with respect to fc in a neighbourhood of 
k = ko for any e, with the possible exception of e = £f; this exceptional instance is (by definition) reserved for Fermi 
liquids (see above). To be explicit, consider ak and f3k which feature in the above expressions for the self-energy of a 
Fermi liquid (see also § V) . Our above statements imply that these may not necessarily be continuously differentiable 
functions of fc in a neighbourhood of k = kp; possible singularities associated with dak/dk and dfik/dk, as fc — > kp, 
are suppressed by {e — ep)"^ and (e — ep) respectively for e = ep. In § VI we establish relationships between the 
behaviour of S(fc; e) for fc and e in neighbourhoods of kp and ep respectively and that of the momentum distribution 
function for fc approaching kp. 

Aside from the above-indicated 'smoothness' assumption with regard to the dispersion of the QP energies which is 
implicit in Luttinger's (1961) proof, there is a second aspect specific to this proof that makes the theorem inapplicable 
to metallic systems of particles interacting through the long-range Coulomb interaction function. Luttinger's approach 
is based on the perturbation expansion of the self-energy operator in terms of the exact single-particle Green function 
(the terms in this expansion are represented by means of skeleton self-energy Feynman diagrams (Luttinger and Ward 
I960)) and the hare electron-electron interaction function. It can be shown (§ IV. D) that a certain class of terms 



non-metallic systems are hy definition non-Fermi liquids. 

* We should like to emphasise that often from the contexts of arguments presented in the literature, in disapproval of the 
many-body perturbation theory, it is not evident what precisely is meant by the 'perturbation theory'. In the present work 
we consider a series for the self-energy operator in terms of the interacting single-particle Green function, each term of which 
is diagrammatically represented by a slceleton diagram, as constituting a perturbation series. This we do entirely in the spirit 
of Luttinger's (1961) work. We should point out that such a series is not a power series in the coupling constant of the 
electron-electron interaction and consequently does not conform to the orthodox definition for perturbation series. 

^ A function f{x) is continuously differentiable in the closed interval [a, b] provided firstly that it is continuous in [a, b] and 
secondly that its derivative exists at all points of the open interval (a, b) and coincides at all such points with a function which 
is continuous in [a, b]; f{x) is continuously differentiable in a neighbourliood of x — xq if xq is interior to a finite interval, such 
as [a,b], in which f{x) is continuously differentiable. 

The notion of a QP without an associated energy dispersion is devoid of physical significance. 
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in this expansion are unbounded as a consequence of non-integrablc singularities corresponding to zero momentum- 
transfer scattering processes (infrared divergence). Therefore, for such systems a formal term-by-term analysis of the 
perturbation series for the self-energy operator can only be meaningful when the perturbation expansion is in terms 
of the dynamically-screened interaction fimction (Hubbard 1957) W{e) (for a comprehensive discussion see Mattuck 
(1992, § 10.4)). In contrast to the static bare Coulomb interaction Vc, W{e) is a function with singularities along 
the real energy (i.e. e) axis (branch points, branch cuts, etc. — similar to the Green function and the self-energy 
operator) and a Lehmann-type spectral representation of this function reveals that these singularities coincide with 
energies of the neutral elementary excitations of the system (i.e. energies of excited iV-electron states as measured 
with respect to the energy of the A''-electron ground state). It follows that an analysis of the self-energy operator along 
the lines of Luttinger (1961) requires additional knowledge with regard to the dispersion of these neutral excitations. 
In other words, in this case the dispersion of the QP energies and that of neutral excitations cannot be considered 
separately, but must be dealt with in conjunction. This is of relevance particularly to uniform two-dimensional systems 
of electrons whose spectrum of coherent neutral excitations (i.e. plasmons) is gap-less, implying some non-negligible 
amount of interference amongst the neutral and the single-particle, i.e. QP, excitations. 

Our present work has been motivated by a long-standing discussion in the literature with regard to the unusual 
properties of layered high- Tc compounds in their normal states. On the one hand these unusual properties have 
been ascribed to the fact that the normal metallic states of these systems not Fermi-liquid states (Anderson 1988, 
1989, Varma, et al., 1989, Anderson 1990a,b, 1991, Littlewood and Varma 1991, Kotliar, et al., 1991, Anderson 
1992, 1993, 1997) while, on the other hand, according to the Luttinger (1961) theorem dealt with in the present 
work, metallic interacting systems in spatial dimensions larger than unity must be Fermi liquids. As emphasised by 
Luttinger (1961), the validity of the Luttinger theorem is conditional to that of the many-body perturbation theory 
for the self-energy operator to all orders. In this light, it is asserted that a non- Fermi-liquid metallic behaviour of a 
system would necessarily imply break-down of the many-body perturbation theory for this system (Anderson 1997). 
These observations give rise to the following questions. 

(1) Can interacting electron systems in spatial dimensions d larger than unity, in particular for d = 2,he non-Fermi 
liquids? 

(2) does the many-body perturbation theory break down in cases where the metallic systems under consideration 

are non-Fermi liquids? 

Concerning question (1), the analyses presented in this work lead us to the conclusion that in principle nothing 
stands in the way of realisation of non-Fermi-liquid states in spatial dimensions d larger than unity. In fact the very 
strict conditions (A) and (B), introduced above, imposed on the self-energy of Fermi liquids, suggests that non- Fermi 
liquid metallic systems (in d > 1) may not be as uncommon as generally perceived. We do not touch upon the above 
question (2) and refrain from statements that are likely to be speculative at this stage. In this connection we point out 
that, since the existence of non-Fermi liquids (in particular for d = 2) and the validity of the many-body perturbation 
theory as applied to these are not a priori mutually exclusive, the question with regard to breakdown of this theory 
is not as urgent as it would be otherwise. This, however, does not diminish the relevance of question (2). Wo should 
like to emphasise that in calculating Y,{k;c), whether perturbatively or otherwise, account has to be taken of the fact 
that any approximation scheme that involves indiscriminate Taylor expansions of functions of energy and momentum 
(which expansions implicitly imply continuous differentiability of these functions in the pertinent regions of energy 
and momentum) may inhibit a non-differentiable E(fc; e) as function of both k and £, and therefore a non-Fermi liquid 
behaviour, from being obtained. 

The above questions, (1) and (2), have been subject of intensive study in recent years, the main body of the results 
thus far obtained pointing towards a Fermi-liquid state in two spatial dimensions. In spite of these, it has as yet not 
proved possible unequivocally to rule out the existence of non- Fermi-liquid states in two-dimensional metallic systems. 
This can be ascribed to two main reasons. 

(i) Owing to the complexity of the many-body problem at hand, application of reliable theoretical techniques must 
of necessity be accompanied by simplifying approximations, the integrity of which may be a matter of dispute. 

(ii) Even relevance of certain calculated quantities to the problem at hand has been matter of debate. 

Both these aspects are aptly represented in the following: Engelbrecht and Randeria (1990) have found no evidence 
for the breakdown of the many-body perturbation theory and the Fermi-liquid theory in a dilute two-dimensional sys- 
tem of fermions interacting through a short-range repulsive potential, contradicting the suggestion made by Anderson 
(1990b). It appears, however, that the phase shift as calculated by Engelbrecht and Randeria is not that which is 
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encountered in the arguments by Anderson (see Engelbrecht and Randeria 1991, Anderson 1991). 

Work by Fujimoto (1990), Fukuyama, Narikiyo and Hasegawa (1991) and Fukuyama, Hasegawa and Narikiyo (1991) 
on the two-dimensional Hubbard model with repulsive on-site interaction U within the i-matrix approximation (which 
is appropriate to the low-density limit) has indicated that, whereas for k ^ kp, ImS(fc,e) retains the conventional 
Fermi-liquid form in three spatial dimensions (presented above), ^ ImI](A; — kp,£)^'{£ — epY \n.\£ — £p\ as e — > £p. 
This singular contribution, obtained earlier by Hodges, Smith and Wilkins (1971) and Bloom (1975) concerning two- 
dimensional Fermi systems, is not sufficiently strong to render the self-energy non-Fermi-liquid like. |^ Serene and 
Hess (1991), using the conserving 'fluctuating-exchange (FLEX) approximation' on a finite but large lattice, have 
also not found evidence for a non-Fermi-liquid behaviour in the two-dimensional Hubbard model. In a more recent 
work, Yokoyama and Fukuyama (1997) have arrived at the conclusion that in the two-dimensional Hubbard model 
with repulsive on-site interaction, the process of forward scattering gives rise to an anomalous contribution to the 
self-energy and consequently to a vanishing quasi-particle weight factor Zkp ■ This result, Yokoyama and Fukuyama 

(1997) asserted, would be "a microscopic demonstration of the claim by Anderson [(1990a, b), (1993)]". We (Farid 
1999a) have shown that the indicated anomalous contribution is entirely a consequence of the violation of a crucial 
symmetry in the momentum space (associated with the time-reversal symmetry of the ground state of the system 
under consideration) by the particle-particle correlation function employed by these workers. 

Momentum-space perturbative renormalisation-group calculations by Shankar (1991, 1994) have also borne out 
a Fermi-liquid picture of fermions in two spatial dimensions. Shankar (1994) enumerated, however, a number of 
possibilities that in principle may render this finding, in its generality, invalid (see § XI in (Shankar 1994)). Work 
by Castellani, Di Castro (1994), Castellani, Di Castro and Metzner (1994), and Metzner, Castellani and Di Castro 

(1998) modelled on the treatment of one-dimensional interacting systems put forward by Dzyaloshinskii and Larkin 
(1973), which exploits the conservation laws and the associated Ward identities, and treats d, the spatial dimension, 
as a real variable, establishes that metallic systems with strong forward scattering (i.e. that corresponding to small 
momentum transfers) are Fermi liquids in d > 1, even though for c? ^ 3 some properties of these systems differ from 
those pertaining to conventional Fermi liquids in d = 3. Here the limited applicability of the formalism to systems 
with dominant forward scattering does not rule out non- Fermi- liquid systems in c? > 1. In a forthcoming work (Farid 
1999b) we present a detailed study concerning some limiting aspects associated with the technique employed in these 
studies. 

Systems of particles interacting through long-range repulsive interaction functions v{r — r') which for large ||r — r'|| 
behave like ^ l/||r — r'jp^'', when 1 < d < 2, and like ~ ln||r — r'||, when d = 2, have been shown to resemble 
one-dimcnsional Luttingcr liquids for 1 < d < 2 and a "Zkp = Fermi liquid" for d = 2 (Bares and Wen 1993, 
Kwon, Houghton and Marston 1995). We note that in metals the Coulomb interaction function, which behaves like 
l/||r — r'll, is screened through the mediation of the charge polarisation fluctuations (Hubbard 1957); within the 
framework of the random-phase approximation (RPA), the screened interaction function in the static limit can be 
shown to behave like ^ cos(2fci?||r — r'||)/||r — r'||'^ (Fetter and Walecka 1971, pp. 178 and 179; Ashcroft and Mermin 
1981, p. 343). It follows that the above long-range interactions should have their origin in processes that are not 
related to the electrostatic electron-electron interaction. The only known interaction that remains long-ranged in 
metals corresponds to that of electronic currents mediated by the exchange of transverse photons (Holstcin, Norton, 
and Pincus 1973, Reizer 1989, 1991). Lack of the static screening of this interaction leads to the asymptotic behaviour 
of the self-energy T,{k; £k) (on-the-mass-shell self-energy), |^ as A; — > fci? (or £k £p), to be that characteristic of the 
marginal Fermi liquids (Varma, et al, 1989, Littlewood and Varma 1991, Kotliar, et al. 1991); on the other hand, for 
the fixed k ^ kp, I](fc;e) behaves Fermi-liquid like as £ ^ £p (Holstein, Norton, and Pincus 1973). We note that the 
S(A:; e) considered here does not involve the effects of the Coulomb interaction beyond a mean-field level; it merely 
describes the self-energy of the current-current interaction whose bare coupling constant is proportional to the square 
of the ratio of the Fermi velocity to the light velocity in vacuum (^ 10^*) and therefore is not of dominating influence 
except in extremely pure metals and at very low temperatures. We shall not enter into further details concerning this 



For a detailed discussion of the singular effective interaction amongst QPs, as encountered in the arguments by Anderson 
(1991), see Stamp (1993); see also Houghton, Kwon, and Marston (1994). 

Throughout, ~' indicates that the corresponding asymptotic relations are correct up to multiphcative constants. 

It can be shown (see Appendix C) that the associated ReE(fcF; e) — S(fcF; eF)~'(£ — £f) so that E(A:f; e) is a continuously 
difFerentiable function of e in a neighbourhood of e = ef and thus, on account of condition (A) discussed above, gives rise to a 
finite Zkp ■ 

For the equation satisfied by the QP energy dispersion £k see Eq. (hsl) below. 
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subject here |^ and refer the reader to the cited Hterature. |^ 

The organisation of this work is as follows. In § II we present the Lehmann spectral representation for the single- 
particle Green function and clarify certain elements that are of particular relevance to our analysis of the Luttinger 
(1961) theorem. In § III we rigorously define the notion of QP and derive equations from which QP energies and 
wavefunctions can be obtained. Here we indicate the necessary steps to be undertaken for obtaining complex-valued 
QP energies. We devote § IV to the main objective of our work, namely a detailed analysis of the Luttinger theorem. 
In § V we briefly review and discuss a theorem due to Migdal (1957) in light of our findings in § IV. In § VI we 
compare Fermi liquids with non- Fermi liquids in terms of the 'smoothness' properties of S(fc; e), as a function of both 
k and e. In § VII we conclude this work by a brief discussion and a review of our main results. In Appendix A we 
introduce and discuss a number of mathematical concepts that we repeatedly encounter in the present work. Here 
we further attempt to clarify the physical relevance of the introduced notions by means of some simple examples. In 
Appendix B we derive the leading asymptotic contribution to the self-energy operator E(e) at large |e|. For this, the 
next-to-the-leading-order asymptotic contribution to the single-particle Green function needs be calculated; in the 
same Appendix we present this contribution as well as some details underlying its derivation. In Appendices C and D 
we consider the next-to-the-leading-order asymptotic term to the self-energy pertaining to Fermi- and marginal Fermi- 
liquids (Appendix C) and the Luttinger liquid (Appendix D). Here we separately deal with the cases corresponding to 
k = kp and k ^ kp. In Appendix D we explicitly demonstrate that the self-energy of the one-dimensional Luttinger 
model is not continuously differentiable at the Fermi points, in full conformity with a general result established in 
Appendix C. In this work we identify electrons with spin-less fermions. 



II. THE SINGLE-PARTICLE GREEN FUNCTION 



Here we explicitly deal with the single-particle Green function. In addition to an exposition of the formal significance 
of the singularities of this fmiction, we discuss three specific and distinct energies fj-, fiN and fiN+i, that are invariably 
(but un-justifiably) identified in the literature concerning the free-electron system. 

Consider the following Lehmann (1954) representation (Fetter and Walecka 1971) for the ('physical') single-particle 
Green function in the coordinate-free representation 



G(s)^?.VA..At| Q(^'"-.' 



irj 



e(£, - ^i) 



e - Es 



-]■ 

irj J 



(1) 



where 



A,:=. 



(^'Ar_i^|V'|*Ar.o), 
(*Ar,o|V'l*W+l.s), 



when Eg < ^, 
when £s > 



(2) 



denotes a 'Lehmann amplitude', and 



En^ — En- 
En+i,s — E 



l,s, 



when Es < 
when Es > /J.. 



(3) 



Above stands for the exact M-particle eigenstate of the interacting system and Em.s for the corresponding 

eigenenergy; s denotes the set of all quantum numbers that uniquely specify |5'a/.s), (we choose s = to signify the 
ground state which we assume to be non-degenerate); ip is the annihilation field operator. |^ Here the 'chemical 
potential' /i is any real value which satisfies < M < l^-N+i, where iJiN'-=EN,o — £^Ar-i,o and iJiN+i'-=EN+ifi ~ En,q. 
That such a /i should exist follows from the requirement of stability of the ground state, meaning that £g:=(i?7v+i,o ^ 
En.o) ^ [En.o ^ £^7V-i.o) = fJ'N+i ^ ^-N be non-negative. For a system with uniform electron distribution, which is 



In particular we do not touch upon the subject of electrons interacting with a gauge field, such as is the case in two- 
dimensional electron systems exposed to external magnetic field (the gauge field here being the statistical Chern-Simons field) , 
where in principle at fractional Landau-level filling factors u with even denominators (specifically at and close to !^ = 1/2) the 
states are metallic but may be non- Fermi liquids (Kalmeyer and Zhang 1992, Halperin, Lee and Read 1993). 
^® For a concise review see (Tsvelik 1995, Ch. 12). 

^"^ Throughout this work, A'^ denotes the actual number of the electrons in the system. 



6 



thus a system in the thermodynamic limit, it holds fiN+i = + 0{N~p), with p > 0. ^ Therefore with regard to 
uniform systems, we are considering a case in which (/xjv, /Ltjv+i) is a finite infinitesimal open interval. In § IV. C we 
shall see that = £f, the Fermi energy of the system of N electrons. 

It can be easily verified that (see Appendix A and in particular § A. 3, for our notational conventions) 

Gizy.^nj:^ (4) 

s 

is the analytic continuation (Whittaker and Watson 1927, pp. 96-98, Titchmarsh 1939, pp. 138-164) of the ('physical') 

single-particle Green function G(e) into the complex z-plane, that it is G(z) on the physical RS (see Appendix A); 
G(e) is obtained from G{z) according to the following prescription 

/(e):=lim7(£±i?7), for e^^. (5) 

For a finite system, the Lehmann representation indicates that G(e) = G{s) has poles |^ at e = for all s (see 
Eq. (^). In the thermodynamic limit, the singular points of G{z) correspond to branch points (isolated singularities 
are not a priori excluded) and continua of branch cuts covering (parts of) the real energy axis (see § A. 4). For finite 
systems as well as those in the thermodynamic limit, G{z) is analytic for all z with lm{z) ^ 0. 

In § III we shall encounter the inverse operator G~^{z). In view of this, we point out that from the representation 
in Eq. (|4|) it can be easily demonstrated that G(z) has no complex zeros (or, better, zero eigenvalues) on the 2;-plane 
(Luttinger 1961). Thus G^^(z), similar to G{z), is analytic for all z with Im(z) ^ 0. 



III. QUASI-PARTICLES; THEIR "ENERGIES" AND WAVEFUNCTIONS 

As we have mentioned in § I, to quasi-particles correspond one-particle wavefunctions which are eigenfunctions of 
some one-particle-like Schrodinger equation corresponding to an energy-dependent non-Hermitian one-body Hamilto- 
nian, the potential-energy part of which consists of the self-energy operator ^ and of possibly a non-trivial contribution 
due to an external potential (such as an ionic potential in a solid). We devote this Section to derivation as well as to 
the discussion of the equations for QP energies and wavefunctions. Some specific properties of the equation for the 
QP energies are of considerable relevance to our discussion of the Luttinger (1961) theorem in § IV. 

In the representation-free form and for the complex energy from the Dyson equation G{z) = Go{z) + 
Go{z)Y,{z)G{z) , with Go{z) the single-particle Green function pertaining to the non-interacting system, one obtains 

G~\z) = l{zl-niz)} (6) 



The value p = 1 as given in, e.g.. Fetter and Walecka (1971, p. 75), is incorrect. 

Below f(z) (./(e)) is a representative for any function of z (e) that we encounter in the present work. 

It should be noted that the prescription given in Eq. (a) is connected to our convention with regard to the time-Fourier 
transform of the time-dependent functions, namely f{e)~J_^dt F{t) exp{+iet/h). Had we chosen exp{—iet/h) rather than 
exp{+iet/h), '±' would have to be 

There is a subtlety involved here. For an open stable system, the spectrum consists of both a discrete and a continuous part. 
The former part corresponds to bound states, and the latter to scattering states; the continuous part of the spectrum gives 
rise to branch-cut discontinuities in the energy-dependent correlation functions of this system. Further, even though in such a 
system the number of particles is finite, the set of discrete one-particle excitation energies has an accumulation point which, 
as indicated in § A.l, is not an isolated singularity of G{z). For closed and finite systems (that is those with impenetrable 
boundaries), the completeness of the one-particle eigenstates implies that the set of one-particle eigenenergies {sa} possesses 
at least one accumulation point (this follows from the Bolzano- Weierstrass theorem — see Whittaker and Watson (1927, pp. 
12 and 13)). Therefore, the non-isolated nature of the singular points of, for example, G{z) is not excJusivefy a peculiarity of 
systems in the thermodynamic limit. 

In spite of the fact that according to our present convention E(fc; e) has dimension s~^, i.e. inverse second, we refer to it 
as self- energy. 
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where / denotes the unit operator in the single-particle Hilbert space and 



n{z):^Ha + nY.{z). (7) 

Any singular point of G{z) is a solution of det (G^^(z)) — 0, or, making use of Eq. (||), of 

det {zl ~ n{z)) = 0. (8) 

This is reminiscent of the equation for the energies of the 'ideal' (i.e. non-interacting) QPs, namely det {zl — i?o) — 0; 
the difference between the two arises from the difference between H{z) and Hq which according to Eq. (|^) is equal to 
h'S{z). The latter, a function of the complex energy parameter z, is a non-Hermitian operator; even for z —>■ s ± irj, 
f] i 0, for £^£_F, leading to I^{z) I](e) (see Eq. (||) above), Hiz) = Hie) is not in general Hermitian. This imphes, 
among others, that, unless 'H{z) is Hermitian, the spectral decomposition of 7i(z) involves the distinct sets of left and 
right eigenvectors, {(pe{z)} and {ipe{z)} respectively. Since the sets of eigenvalues pertaining to the latter two sets 
are up to ordering identical (Morse and Feshbach 1953, p. 885), by choosing the appropriate ordering, the spectral 
(or bi-orthonormal (Morse and Feshbach 1953, pp. 884-886, Layzer 1963) — see later) representation of T-C{z) is as 
follows 

H(z)=^^,(z)V^,(z)0l(z), (9) 

where Ei(z) denotes the common eigenvalue corresponding to (j)e{z) and ijjeiz)'. 

i\{z)n{z) = Et{z)i\{z)- U{z)^i{z) = E,{z)Mz)- (10) 

For Et{z) ^ Ei'{z) it holds {(pi{z),ipe'{z)) — Si^t (Morse and Feshbach 1953, p. 885). In the case of degeneracy, 
that is Et{z) — E('{z) for £ ^ the degenerate left and right eigenvectors can be made orthogonal through a 
Gram-Schmidt orthogonalisation procedure. When the left and right eigenfunctions are simultaneously bases of the 
unitary irreducible representations of the symmetry group of the QP Schrodinger equation, Eq. (|l^) , the degenerate 
eigenvectors are automatically orthogonal^ as long as _they belong to different unitary irreducible representations 
(Cornwell 1984, pp. 81-83). |^ Further, {(l)i{z)} and {'ipi{z)} satisfy the completeness relation (Morse and Feshbach 
1953, p. 886) 

J2Mz)4iz) = i, (11) 

e 

or V^f (r;^)<^|(r'; z) = (5(r - r'). 

Because H'^iz) = H{z*), for Im(z) 7^ 0, a property that follows from S^(z) = ^{z*) (Appendix B in DuBois (1959a), 
Luttinger 1961), it can be readily shown that, for Im(z) 7^ 0, 

Mz) = Mz*), Mz)=Mz*), Et{z*) = E}{z). (12) 

As we shall see later, the last result ^ is of particular significance to our considerations. By defining Ei{£) in 
accordance with Eq. (||), for Imi?^ (e) the following must hold 



Consider, for instance, the uniform-electron system with which we deal in detail in the present work. The symmetry 
group of this system is the continuous translation group (which is Abelian) and therefore the corresponding bases for the (one- 
dimensional) unitary irreducible representations are specified by wave-vectors k. The time-reversal symmetry of the problem 
implies degeneracy of eigenfunctions at k and — k ("Kramers' degeneracy" (Landau and Lifshitz 1977, pp. 223-226)). Since 
for k 7^ 0, k 7^ — k, the eigenfunctions corresponding to k and — k are therefore automatically orthogonal. Note that, since 
our 'electrons' are spin-less, for the application of the Kramers theorem we cannot rely on the condition that the total spin 
of the system is half-integer. However, we can rely on the fact that all irreducible representations of the translation group 
corresponding to k 7^ are essentially complex. 

The left-hand side of this relation is the so-called idemfactor of the one-particle Hilbert space. 

This result is the expression of the Riemann-Schwarz reflection principle (Titchmarsh 1939, p. 155) which follows from the 
analyticity of ^(2) over the entire complex z-plane, with the exception of the real e-axis, and the fact that over the finite (even 
though infinitesimal) open interval (^jv,pjv+i), S(e), or G{e), is real valued. 



8 



ImEiie) I 0, for e > ^. (13) 

Violation of these inequalities signifies breakdown of causality, or instability of the ground statejiue to its collapse into 
a lower-energy state; to appreciate this, note that (combine Eqs. (^) and ( [Til) ) G{z) — h,^^%})i{z)(t)\{z) / {^z — Ei{z)^ 
and compare this with the Lehmann representation in Eq. (|l|) (see also Eq. (|^)). 

From the spectral representation for H^z) in Eq. (^ above, together with Eq. (11), the equivalence of Eq. (||) with 
the set of equations 

E,{z)^z, ye, (14) 

is readily established. As we have mentioned in § II, singular points of G(e), according to the Lehmann representation, 
coincide with the real- valued quantities £s (see Eq. (^). Since Eq. ( p^ is the equation for these singular points, it 
may be expected that by taking the appropriate limits z — > eizyy, with ij I (see Eq. (^)), Eq. (14) should transform 



into an equation with real-valued solutions e^; owing to the analyticity of G~^{z) for Im(z) ^ (see § II), it follows 
that Eq. ( p^ ) cannot be satisfied for any z with Im(z) ^ (for a physical interpretation of this statement see the 
following paragraph). As we shall demonstrate below, in the thermodynamic limit the number of real- valued solutions 
of Eq. (p^), with z — > e ± 177 (sec Eq. (||)), is limited. Further, in the same limit, these real-valued solutions are in 
general not poles, in contrast with what the Lehmann representation in Eq. (^ (or Eq. (^)) may suggest. This is not 
a contradiction, since as we indicate in § A. 4, analytic properties of correlation functions in the thermodynamic limit 
are fundamentally different from those of finite systems. 

It is often mentioned that finite lifetimes of QPs is reflected in the energies of these being complex valued. This 
statement, without further qualification, is misleading. To clarify this, consider Eq. (Wm and suppose that it were 



satisfied by z = zq with Im(zo) 7^ 0. Owing to the reflection property of Ei{z) presented in Eq. (12), it follows that 
z = Zq must also satisfy Eq. (p^). This cannot be the case, since depending on whether Re(zo) is less than or larger 
than /i, one of the two solutions zq and Zq signals violation of the principle of causality. In this connection, recall that 
the single-particle Green function has been defined in terms of the time-ordered product of two field operators in the 
Heisenberg representation (Fetter and Walecka 1971). Consequently, no zq can satisfy Eq. ( p^ ) unless Im(zo) = 0. 
Thus for a given £, Eq. (^4|) has either a real-valued solution or no solution at all; the possible real-valued solutions 
(see § A.l) are, in general, not isolated. To obtain complex-valued quasi-particle energies, the appropriate equation 
to be solved is the following 

Ee{z) = z, V^, (15) 

where Ee{z) represent the unique analytic function of which Ei(z) is one specific branch (see § A. 3). The analytic 

function Ee{z) (i.e. any of its branches; there may be infinity of these, depending on the nature of the branch points of 
Ei{z)) can in principle be calculated from the knowledge of Ei{z). To this end, one starts from E£{z) and analytically 
continues it across its branch cuts into a RS 'above' or 'below' the physical RS (see examples in §§ A. 3 and A. 4). This 

process can be repeated on any of the RSs and in this way one recovers all the possible branches of E{z). From the 
specification presented in Eq. (^, it is readily observed that, to obtain the physically most relevant complex- valued 
solutions of Eq. (^5|), one has to analytically continue Ei{z) from the first (third) quadrant of the physical RS into the 

fourth (second) quadrant of a non-physical RS when Re(z) > /x (Re(z) < /z). We note that solution zq of Ei{z) = z 
with Re(zo) > ^ corresponds to a QP that propagates and attenuates forwards in time (i.e. it concerns a 'particle-like' 
QP) and that with Re(zo) < /i to one that propagates and attenuates backwards in time (a 'hole-like' QP). 



IV. DISCUSSION OF THE LUTTINGER THEOREM 



We are now in a position to deal with the central objective of the present work. In doing so we distinguish three 
major steps in the proof by Luttinger (1961). We shall separately describe each of these steps and, when necessary, 
comment on the details. We conclude that Luttinger's theorem does not exclude the existence of non-Fermi-liquid 
metallic states for systems in spatial dimensions larger than unity. Rather, it merely establishes that Landau's 
Fermi- liquid theory for metals is consistent with the principles of the many-body theory of interacting systems. 
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A. Preliminaries 



Before proceeding with the proof of the Luttinger (1961) theorem, several remarks are in order. 

First, our considerations in § II have made evident that /ijv and fiN+i sue two branch points of G{z) and 2(2) and 
consequently of Ei{z), for all (.. For simplicity of notation, let f{z) denote any of these functions of z. The fact that 
UN and fJ-N+i are branch points of f{z) implies that such expansion as lmS{k;e) ~ =Fafc(e — Ef)^, e^ef, cannot 
correspond to the leading-order term in a Taylor expansion of S(fc; z) when z — ep identified with either /xjv or 
fiN+i- Rather, it corresponds to the leading-order term in an asymptotic expansion (Whittaker and Watson 1927, 
Ch. VIII, Copson 1965, Dingle 1973, Lauwerier 1977) of ImS(fc; e) for e ^ (see § A.2). 

In view of the fact that for the uniform-electron system under consideration, the interval {fJ-N, I^n+i) is infinitesimally 
small, it is tempting to identify ep, fJ-N and fJ-N+i, as is commonly done in the literature. Such an identification is not 
justified. Our discussion of a Migdal's (1957) theorem in § V further clarifies this statement. We note that, strictly 
speaking, ef is identical with /ijv- Thus the precise statement of the Luttinger theorem under consideration is as 
follows. ImE(A:;£) = for e e {^jy, fiN+i)', ImS(fc;e + irf) = Iml](fc;£) ~ — Q!fe(e — ^.n+iY fo'' ^ ^ MJV+i and 
ImE(fc; £ — irj) = Iml](fc; e) ^ ak{s — ^jln)"^ for ^ < I^n- 

Second, the above asymptotic expansions are not uniform (see § A.2), as is evident from the sign of ImS(fc; z) which 
depends on the sign of Im(z) (or on arg(z)). Thus z — fj,N and z — ij-n+i are indeed branch points. 

Third, analysis of the next-to-leading order terms in the asymptotic expansions of I](fc; z) around z = fj.^ and 
z — fJ-N+i reveal that these terms involve the logarithm function, implying that both fi^ and fiN+i are branch 

points of infinite order, that is to say f{z) (see § A. 3) consists of infinitely many branches. Already the fact that /ijv 
and /ijv+i are not regular singularities of f{z) suggests that care should be exercised in defining the QP weights, in 
particular Zk^. which features in the Migdal (1957) theorem and which determines the value of the discontinuity in 
the momentum distribution functions at the Fermi momentum hkp. We shall return to this point in § V. 

In spite of our first remark above, for conformity with the existing texts in what follows we shall use the conventional 
notation (such as that in our Abstract) and proceed with the demonstration of the implicit assumption involved in 
the proof of the Luttinger (1961) theorem. To this end, we first summarise the main aspects of Luttinger's (1961) 
proof. This proof is based on the assumption of validity of the many-body perturbation theory and further rests on 
the analyses of the terms in an implicit perturbation expansion of the self-energy operator; the implicit nature of 
this expansion is associated with the fact that contributions in this are in terms of the single-particle Green function 
pertaining to the interacting system which is, in turn, a functional of the self-energy operator itself. Diagrammatically, 
the terms in this perturbation series are represented by the proper skeleton self-energy diagrams (Luttinger and Ward 
1960), that is those connected self-energy diagrams that do not accommodate any simpler self-energy parts in them. 

In the proof by Luttinger (1961), the skeleton diagrams are in terms of the isotropic bare electron-electron interaction 
function v{r — r') and this, as we shall discuss in § IV. D, implies that Luttinger's approach is not directly applicable 
to cases where v(r — r') is the long-range Coulomb interaction function. 

Since the many-body perturbation theory underlies the proof by Luttinger (1961), the Luttinger theorem is valid 
insofar as this perturbation theory is valid. This aspect leads to the natural conclusion that metallic non-Fermi-liquid 
systems in spatial dimensions d larger than unity (see footnote ^ further on) fall outside the domain of applicability 
of the many-body perturbation theory, or that this theory should break down when applied to such systems. As has 
been emphasised earlier (Anderson 1993), validity or otherwise of a many-body perturbation series cannot be decided 
solely on the basis of its convergence or divergence, respectively. Despite convergence, the calculated limit may be 
unphysical. Q 



In the work by Luttinger (1961), there is no explicit mention of the value or range of values of k for which ImE(fc; e) ~ 
=Fafc(e — ef)^, e^Ef, is valid. This is, however, implicit in Eqs. (59)-(61) of Luttinger's work: Eq. (61) implies u to be small in 
magnitude {u:={jj, — x); in our notation, u:—{ef — e)), so that by Eq. (60) one deduces ti, t2 and to be small, requiring by 
Eq. (59) that fci, fe, fca ~ kp- From ka = ki + k2 — k, using a geometrical construction, it is easily deduced that < k^Skp- 
^"^ The following example due to Simon (1970) should be clarifying: Consider the hydrogen-like Hamiltonian in three dimensions 
(here we use the Hartree atomic units): H:— — — A/r, r > 0. The discrete energies of this Hamiltonian for A > are 

£n(A) = — A^/[2n^], n = 1, 2, . . .. For A < there are no discrete levels (no bound states). However, £n{X) being a finite-order 
polynomial (technically, an entire function) of A, the Rayleigh-Schrodinger perturbation expansion for, say, the ground-state 
energy (i.e. £'i(A)), with —A/r playing the role of the 'perturbation', yields exactly — A'^/2, irrespective of whether A > or 
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B. Explicit analysis of the Luttinger theorem 



After indicating that the perturbation contributions in the above-indicated expansion for E(fc;£) in terms of the 
skeleton self-energy diagrams below second order in the electron-electron interaction v are e-independent, Luttinger 
(1961) considered in detail the second-order self-energy diagram that we have depicted in Fig. 1 (for facilitating direct 
comparison with the results in the work by Luttinger, below we adopt the notation that he employed as closely as 
possible). In doing so, Luttinger first evaluated the imaginary part of the self-energy operator corresponding to this 
diagram, with the single-particle Green function pertaining to the interacting system, that is G{e) (represented by 
full line in Fig. 1), replaced by that pertaining to the non-interacting system, i.e. Go(e). ^ From the well-known 
rules for transcribing Feynman diagrams into mathematical expressions in the momentum representation (Fetter and 
Walecka 1971, pp. 100-105), it is readily seen that the expression corresponding to the diagram in Fig. 1 involves three 
wave- vector integrations, over ki, k2 and ka. Following the fact that for the matrix elements of the electron-electron 
interaction, (k, k3|i>|ki, k2) = v(k — ki)^k-ki,k2-k3 holds, where k stands for the external wave- vector (i.e. the wave- 
vector associated with fc:=||k|j in the argument of Ti{k] s)), it would be tempting to eliminate one of the wave-vector 
integrals. This elimination would make the subsequent algebra extremely cumbersome. Rather, Luttinger (1961) 
transforms the three wave- vector integrals into three one-dimensional integrals over energies , and e^^ of the 
non-interacting electrons, where 

el:=^k\ (16) 

with rrif, the bare-electron mass. |^ 

A most crucial ingredient in Luttinger's proof consists of the recognition that, although the momentum conservation, 
enforced by the above ^k-ki,k2-k3, uniquely determines one of the three wave- vectors in terms of the other two, say 
ka in terms of ki and k2, nonetheless, barring one-dimensional systems (Luttinger 1961, footnote 5), e^_^, e^^ and e^^ 
can be independently and continuously varied over some finite range of values. |^ This can be seen as follows: using 
the above dispersion for e° as well as ka = ki + k2 — k, it is readily seen that 



-0 _ , , 
-/C3 — + ^ki "T ^fca 



— {ki •k2-k-ki -k.k2}. (17) 



For fixed values of e^.^ and e^^, and in spite of ka = ki+k2— k, for systems in spatial dimensions d > 1, the orientational 
freedom of vectors ki and k2 with respect to each other as well as k allows for continuous finite variations in e^^ . 
A crucial element in our arguments that follow, is derived from the fact that freedom in the independent variations 
of is a consequence of the expression in Eq. (|lj), which, in turn, follows from the speciGc form of the energy 
dispersion pertaining to non-interacting electrons in Eq. (p^). As a matter of course the latter form is not unique in 
that one can construct energy dispersions different from that in Eq. ( p^ ) that equally give rise to the possibility 
of independent and continuous variations of e^., j = 1,2,3, over some non- vanishing region. At the same time, it is 



A < and independent of the magnitude of |A|; the coefficients of A™', for m = 0, 1 and m > 2, in this expansion are identically 
vanishing. Evidently, the convergence of this perturbation series for A < is a 'bogus' convergence. 

Explicitly, we have Go(fc; e) = h{e{kF - k)/{e - - ir?) 4- e(fc - fcF)/(e - ?7 i 0; el is defined in Eq. (^. 

The approach by Luttinger (1961) aims at obtaining the leading-order term in (e — ef) (in Luttinger's notation, u), so 
that other integrals which merely determine numerical pre-factors are not explicitly dealt with by him. For clarity, consider 
a three-dimensional non-interacting system with the one-particle energy dispersion presented in Eq. (^) . In the spherical- 
polar coordinate system {k,(f),6) we have d^k — k^ sm{6) dk d(f) AO = ^{2me/h'^)^^'^ y^e'lsm{9) de1d(l>d9. In the treatment by 
Luttinger (1961), d^fc cx de^ plays an all-important role. 

This is the instance where one-dimensional interacting systems are singled out as being in some fundamental way different 
from higher-dimensional interacting systems. The Luttinger model, which turns out to be akin to a large class of one-dimensional 
models (Haldane 1980, 1981), is however introduced in 1963 by Luttinger (1963) and was first correctly treated by Mattis and 
Lieb (1965). An earlier model to which the Luttinger model is similar, is due to Tomonaga (1950). For three different treatments 
of the Tomonaga-Luttinger models see Dzyaloshinskii and Larkin (1973), Luther and Peschel (1974), and Everts and Schulz 
(1974). For reviews concerning one-dimensional systems of interacting fermions see Solyom (1979), Mahan (1981, § 4.4), Voit 
(1994), Schonhammer (1997) and Schulz, Cuniberti and Fieri (1998). 
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not difficult to put forward energy dispersions that do not allow for such an independent continuous variation over a 
finite domain. We shall return to this point in § IV. C. 

It is from the observation with regard to the possibility of independent variations of the three energies s^., j = 1,2,3, 
that Luttinger (1961) has established the leading-order contribution to ImE(fc; e) of the second-order diagram in Fig. 1 
to be proportional to {e ~ Sp)'^ as e — > e^. Q 

The second step and, from the viewpoint of our present work, the most crucial step in Luttinger's (1961) proof 
consists in establishing that, upon employing G(e) rather than the G'o(e) of the first step, the leading-order contribution 
to lmS{k; e) pertaining to the second-order self-energy diagram in Fig. 1 is also proportional to (e — ep)'^, for e ^ ep- 
This conclusion is of fundamental importance to the third part in Luttinger's proof, where, after having established 
the latter property, namely that the leading-order contribution to ImE(fe;e), as e ^ e^, is unaffected by evaluating 
the contributions of the skeleton self-energy diagrams in terms of Go{e) rather than G(e), Luttinger demonstrated that 
to all orders in the perturbation theory ImS(fc;e)~'(£ — sp)^ or explicitly, that any skeleton self-energy diagram (in 
terms of Go(e)) of second and higher order in the bare electron-electron interaction v has a contribution to Iml](fc; e) 
proportional to (e — e^.)^™, where m > 1, for e ^ e^i^. ^ 

In this second part of his proof, Luttinger (1961) proceeded as follows. On evaluating the contribution of the 
diagram in Fig. 1 to ImS(A;;e) in terms of G{e), he made a simplifying approximation, which, roughly speaking, 
corresponds to replacing a Lorentzian by a Dirac (5- function. ^ Upon this, Luttinger arrived at an expression which 
has the same formal structure as the one discussed above. Similar to the second-order expression for ImS(fc; e) in terms 
of Go, the simplified expression in terms of G involves three wave- vector integrals over ki, k2 and k^. Now, at this 
place, Luttinger took a step which, as will become evident below, implicitly contains the property to be proven. Here 
Luttinger transforms the three wave- vector integrals into three energy integrals, the energy dispersions being Sk-i, Sk2 
and £fe3 , to be compared with , and considered above (see § IV. C). Note that the dispersions Sk^ , = 1, 2, 3, 
correspond to the interacting (i.e. true) QPs and satisfy Eq. ( |l^ or Eq. (^sj) below. According to Luttinger, 
all the necessary mathematical steps that are to be taken in order to obtain, from the last-indicated integral, the 
leading-order contribution to ImE(fc; e)^ as e ^ sp, are identical with those taken in dealing with ImE(A:; e) evaluated 
in terms of Gq. This is an unjustified statement. The dependence on k of £/c being unknown, there is no a priori 
reason why can be varied continuously and independently of and Efcj over a finite range. If the effect of 
electron-electron interaction on the dispersion of the QP energies close to the Fermi surface were only to re-normalise 
the QP mass, then of course for A; in a neighbourhood oi k = kp one invariably had Sk ?i^fc^/(2m*), with m* an 
electron's renormalised mass (Pines and Nozieres 1966, Abrikosov, Gorkov and Dzyaloshinski 1975). In such an event, 
for fc in a neighbourhood of k = kp, Eq. ( p^ ) would hold for with me replaced by to*, and thus the just-mentioned 
assertion by Luttinger (1961) were correct. However, as we shall demonstrate in § IV. C, Sk ~ ?i^fc^/(2TO*) can only 
apply if S(fc;e) is of the Fermi-liquid type, thus establishing that Luttinger's proof amounts to a demonstration of 
consistency of the Fermi-liquid state with the many-body theory of interacting metals. This demonstration does not 



We note in passing that 'it^' in Eq. (60) of Luttinger (1960) is in error; the correct pre-factor is unity. 

It is known (Hodges, Smith and Wilkins 1971, Bloom 1975, Fujimoto 1990, Fukuyama, Narikiyo and Hasegawa 1991, 
Fukuyama, Hasegawa and Narikiyo 1991) that, in isotropic two-dimensional interacting systems, ImE(fc = kF;e)^'{e — 
£f)^ In |e — £f| as e ^ ef- Since the leading e-dependent contribution to ReE(fc = /cf; e) associated with this ImE(fc = fcf; e) 
is proportional to {e — ef) (see Appendix C), we have that E(fcF;e) is a continuously difFerentiable function of e in a neigh- 
bourhood of e = Ef- In view of our statements in § I (see also § IV. C), the above logarithmic contribution does not turn the 
system into a non- Fermi liquid. 

This approximation neglects the finite lifetimes of the QPs. From ImE(fc; ef) = (note in passing that this together with 
e1 + hT,{k;EF) = Ef define the Fermi surface — see footnote ^ below) and the requirement of continuity of ImE(fc; e) for 
E Ef, we have ImE(fc; e) — > for e ef- Consequently, for e — > ef life-time effects are not of relevance to the ieading-order 
contribution to ImE(fc; e) as e ^ ef- The mentioned 'Lorentzian' has its origin in the finite lifetimes of the QPs. For a summary 
of various ways in which ImE(fc; e) can vanish as e — > ef, see the text following Eq. dsi] ) in § V. 

In Luttinger's (1961) work there is no symbolic distinction between, for example, what we have denoted as Y:{k;z) and 
E(fc; 2). Further, Luttinger does not discuss whether Eq. jl^ ) can be exactly satisfied or not (for details see § IV. C). 

The contribution of the diagram in Fig. 1 to ImE(fc;e) for one- dimensional metals, according to Luttinger (1961, p. 946, 
footnote 5) amounts to ~'(e — £f), which we know to be characteristic of marginal- Fermi liquids (see Appendix C) rather than 
one- dimensional Luttinger liquids (see Appendix D) over a physically-relevant range of values for the anomalous dimension 
a = 270, that is < a < 1. At the root of this incorrect inference by Luttinger (1961) lies the fact that in one- dimensional 
metals, Ek, contrary to Luttinger's implicit assumption, is not a continuously-differentiable function of k and therefore d^fc oc de^ 
does not apply, while d^fc oc de^ is valid. See § IV. C. 



12 



rule out existence of states different from the Fermi-liquid state. 



C. On the dispersion of the QP energies 



Here we analyse some general features of the QP energies pertinent to a uniform-electron system. 
Consider Eq. ( p^ specialised to a uniform system. Any possible reaJ- valued QP energy Ek is a solution (which for 
an arbitrary k may not exist; see § V) of the following equation (see also Eq. (26) below) 



(18) 



The fact that E(fc; z) has a branch point at z 
satisfied by the reaJ- valued energy ep (below £%<:—£ 



-£kp (§ IV. A), implies that the following equation is well 



£f = £p + h'E^kp; ef)- 



(19) 



Assuming S(fc;£) to be a continuously differentiable function of k and e in some neighbourhoods of k — kp and 
e — ep, respectively, from Eq. (|l8|) together with the manifest differentiability of e° in Eq. ( p^ ) with respect to k it 
follows that Ek is similarly a continuously differentiable function of A: in a neighbourhood of k — kp. Consequently, ^ 



£fe = ep + hvp{k — kp) + o(^k ~ kp) , 



where 



vp:=- 
n 



1 9£fe 

dk 



k—kfi 



Tikf 



ml 



stands for the Fermi velocity. From Eq. (|l8|) it can be shown that 



dsk 
dk 



^Zk 



d{el + h^{k;ep)} 



k—kf 



where 



Zk,:=[l-h 



dk 



d^{kp;e) 



k—kp 



de 



(20) 



(21) 



(22) 



(23) 



is, according to a Migdal's (1957) theorem (see § V), exactly the amount of discontinuity in the momentum distribution 
function n{p) at the Fermi momentum pp = fikp. Eqs. (|2^) and ( |2l| ) are readily seen to imply that for k ^ kp, 
£k ~ h?k^/(2ml). The above steps make explicit how the latter dispersion relation is rooted in the assumption 
of continuous differentiability of I](fc;efc) in a neighbourhood of fc = kp. For Fermi liquids < Zkp < 1 holds, 
where Zkp = 1 corresponds to the case of non-interacting QPs. The condition Zkp = 0, which is a signature 
of marginal- and Luttinger-liquid systems, corresponds, according to Eq. (p3|), to the case where d'S{kp; e)/de is 
unbounded at e = ep. In cases where Zkp is vanishing, Eq. (^ ) formally implies that, for the Fermi velocity vp 
to be finite, dT,{k;ep)/dk is also to be unbounded &t k = kp. Although formally Eq. (|20| ) can be maintained with 
Vp = —h~^dYi{k;eF)Jdk\k=kp/dYi{kp]e)/ds\e=ep^ it should be evident, however, that in the case at hand the right- 
hand side of Eq. ( po|) neglects the fact that E(fc;e) is singular at both k = kp and e — ep with the singularities 
showing up already in the first derivatives, thus invalidating this formal definition of the Fermi velocity. To appreciate 
the significance of this observation, note that the key element in obtaining Eq. ( po[ ) has been the assumption of 
continuous differentiability of T,{k;e) in some neighbourhoods of k = kp and e — ep, whereby it has been possible 



Isotropy of the system under consideration implies kp to he independent of the electron-electron interaction. We note in 
passing tliat this is a corollary to another celebrated theorem due to Luttinger (1960); see also Luttinger and Ward (1960). 
^'^ As indicated earlier in this work, o in f(x) = o(^g{x)) signifies f{x)/g{x) for x ^ 0. 

Recall that < Zkp < 1 although necessary, is not sufficient for rendering a metallic system a Fermi liquid (see § I, 
conditions (A) and (B)). 
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to write E(fc; sj?) ~ 'E{kF;SF) + P{k — kp) + o{k — kp) and ^{kp; e) ~ Y:{kp;sp) + ^{e — ep) + o(e — ep), with 
/3 and 7 bounded, in some neighbourhoods of k = kp and e — ep, respectively. These asymptotic relations signify 
the facts that {S(fc; sp) — T,{kp; sp)}/ (fc — kp) and {Ti(kp; e) — Yi{kp] ep)} / (e — ep) are vanishing for k kp and 
e ^ respectively. With, for instance, T,{k;ep) ^ T,{kp;eF) + P\k — kp]"^", < 70 < 1 (in the present case 
a:—2"f() is referred to as the anomalous dimension; see Appendix D, in particular paragraph following Eq. (D30)), 
for k kp, or 'E{kp;e) ^ 'S{kp;ep) + 7(2 — eF)ln|e — ep\, for e ep, which have not been accounted for by 
Eq. (|20|), {l^{k;ep) — l^{kp;ep)}/{k — kp) and {T,{kp;e) — ^{kp; ep)}/ {e — ep) clearly diverge for k kp and 
e ep, respectively, and consequently both Eq. ( |20| ) and the notion of Fermi velocity of QPs in the vicinity of the 
Fermi surface become altogether meaningless. 

The following two observations should be clarifying: first, any positive non-integer 70 in the above expression would 
signal a non-analytic singularity in the fc-dependence of S(/c;e), even when 70 > 1; second, even if one puts aside 
the possibility of non-differentiability with respect to k of 'E{k;e) in a neighbourhood of fc = kp, it is certain that 
'E{k;z) possesses an essential singularity at z = = ^i^. For instance, Eq. ( pO| ) with vp defined in accordance 
with Eqs. ( |2l| ) and (22) takes no account of the possibility of, for example, Im5](A;; e)^'{e — ep) whose corresponding 
ReS(A;; e) — l^{k;ep)^'{e — ep) ln\e — ep\ as e ^ ep (see text following Eq. ( ^l| ) below as well as Appendix C) which 
is specific to 'marginal' Fermi liquids (Varma, et al. 1989, Littlewood and Varma 1991, Kotliar, et al. 1991). 

Having demonstrated the difficulties that arise as a consequence of the fact that E(/c;e) (and, owing to Eq. (p^), 
ek) is not a priori 'smooth' (in our above-indicated sense), it becomes evident that the implicit assumption in the 
proof of the Luttinger (1961) theorem, as though T,{k; e) were invariably 'smooth', has no theoretical justification. In 
particular, we should like to emphasise that our above exposition makes evident that no theoretical treatments that 
aim to expose the behaviour of E(A:; e), for e ep, whether S(fc; e) is evaluated by perturbative or non-perturbative 
techniques, should rely on the expansion in Eq. (^), since any consistent treatment that relies on this expansion is 
bound to arriving at the conclusion that S(fc; e) is Fermi-liquid like (validity of the expansion in Eq. ( po[ ) is necessary 
and sufficient for the satisfaction of conditions (A) and (B) in § I); a different conclusion must of necessity signal some 
inconsistency (or inconsistencies) in the treatment, including plain algebraic errors. 

One of the main conclusions that may be drawn from our above analyses is that a non-Fermi-liquid-type S(/c; e) 
does not necessarily signal breakdown of the many-body perturbation theory. Consequently, we observe that non- 
Fermi-liquid behaviour can in principle be manifested in any spatial dimension higher than d = 1. 



D. A specific feature 



Here we comment on one of the specific features in the treatment by Luttinger (1961) which brings out some 
additional aspects concerning the above-indicated implicit assumption. The diagram in Fig. 1 has the property that 
it involves a so-called 'polarisation part', that is a part which can be disconnected from the self-energy diagram 
through removing two interaction lines (indicated by broken lines). Momentum conservation implies that to both of 
these interaction lines the same momentum-transfer vector is associated. Consequently, any polarisation part in a 
self-energy diagram contributes at least a square of the electron-electron interaction to the respective integrand. As a 
consequence, the contribution of the diagram in Fig. 1 is divergent when the electron-electron interaction is the bare 
Coulomb interaction (for a comprehensive treatment see Mattuck (1992, § 10.4)) for which ^(r — r') oc l/|lr — r'|| holds 
and thus v{k) oc 1/fc^ in three spatial dimensions. Thus Luttinger's analysis is not directly applicable to systems 



Compare with Eq. (^o|) . See also footnotes ^ and ^ 

We note in passing that for metallic systems the locus of the k-points satisfying + ^S(^i £f) = £f (compare with Eq. (p^) 
is the Fermi surface, with e2 + ?!-S(fc;£F) < £f (> £f) defining the interior (exterior) of the Fermi sea (Galitskii and Migdal 
1958, Luttinger 1960, Eqs. (6) and (94)). Clearly, the existence of a Fermi surface is independent of the value of Zkp, as the 
latter is determined by the derivative with respect to e of E(/ci?;e); see Eq. (^3|). In particular Zi^^ = does not rule out a 
Fermi surface. We point out that ImE(A;;£F) = 0, for all k (see footnote ^ above), and that in an isotropic system the Fermi 
surface can consist of disconnected concentric surfaces; since is a monotonicalty increasing function of k, such Fermi surfaces 
can exist (in isotropic systems) only in consequence of the electron-electron interaction: for -\- TiT,{k;eF) = £f to have more 
than one solution, it is necessary that 7iE(fc;eF) can counter ej?. In the present work we explicitly deal with a single Fermi 
surface of radius kp- 
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of particles interacting via the Coulomb interaction. However, as is well known, the contribution of the set of 
all divergent polarisation diagrams (the random-phase approximation, RPA, bubble-like diagrams) can be exactly 
calculated (as these diagrams give rise to a geometric series) and one observes that screening effects in the static limit 
remove the singularity of the bare Coulomb interaction function v{k) at k — (see, e.g., Fetter and Walecka 1971, 
pp. 178 and 179). In analysing the leading-order term in ImI](A;;e) as e ^ £f, the e-dependence of the dynamically- 
screened interaction function (Hubbard 1957) W{e) necessitates knowledge not only of the energy dispersion of the 
QP excitations, i.e. e^, but also of that of the bosonic neutral excitations, i.e. e/c; here for the energies of the 
'neutral' excitations we have es-=EN.s — En^q > (Fetter and Walecka 1971). For an interacting system, the exact 
dispersion of is unknown, similar to that of e^; e^ can be in principle determined from the dynamical density-density 
correlation function x(fc; e); for et this function plays a similar role as G(fc; e) does for Sk- It follows that for metallic 
Coulomb systems, the behaviour of ImE(fc; e) cannot be solely determined from the perturbation series expansion for 
the self-energy operator; one needs in addition the perturbation series expansion for x{k;e). 

From the above arguments we conclude that, even if the implicit assumption by Luttinger (1961), discussed in 
§§ IV. B and IV. C, were correct, Luttinger's asymptotic expression for ImE(fc;e) as presented in the Abstract, would 
not be justifiably applicable to self-energies of metallic systems of particles interacting via the Coulomb interaction 
function (sec § I for references to studies on systems of particles interacting through long-range repulsive interaction 
functions). 

V. A BRIEF DISCUSSION OF A MIGDAL'S THEOREM 

In § IV we have mentioned that the quantity Z^p , which coincides with weight of the Landau quasi-particles on the 
Fermi surface, is also equal to the amount of discontinuity of the momentum-distribution function n{p) at pp = fikp. 
This statement constitutes a theorem due to Migdal (1957) (Luttinger 1960). Here we briefly discuss this theorem 
in light of our above considerations. In what follows we deal with n(fc) and use the commonly-employed designation 
'momentum-distribution function', despite the fact that k =p/'h is not momentum. For n(fc) we have 

1 /■ dz ~ 

y^{k)■.={^N.o\K^k\'^N.o) = ^ G{k-z), (24) 

where a\. and au stand for creation and annihilation operators, respectively. The contour C of integration is depicted 
in Fig. 2. To prevent confusion, we point out that the right-hand side of Eq. ( p^ differs by a minus sign from the 
corresponding expression in Migdal's (1957) work. This is because the single-particle Green function adopted here is 
defined according to the modern convention and is equal to minus the Green function employed by Migdal (1957). 
Further, our use of the symbol G{k; z) (rather than G{k] e) which is the common notation) is in accordance with our 
present conventions (see § A. 3) and has the advantage that our following discussions will not suffer from mathematical 
ambiguities. One of these ambiguities can be found in the original work by Migdal (1957). Migdal stated namely 
that, for k infinitesimally less than kp (denoted by A:]^), G{k; z) would have a pole, with an infinitesimal imaginary 
part (apparently due to ImS(fc;e) as indicated in the Abstract), enclosed by C. Migdal (1957) further asserted that, 
through increasing k from kp to kp {kp denotes a value infinitesimally larger than kp), the imaginary part of the 
mentioned pole would change sign, upon which this pole would leave the interior of C. Since the regular or incoherent 
part of G(k; z) does not contribute to the difference n(fcj^) — n{k'^), from Eq. ( p^ ) it would follow that Z^^ would be 
the residue of an isolated pole of G{kp; z) at z — ep. The Migdal theorem (Migdal 1957, Luttinger 1960) thus states 

nikp)-n{k+)^ Zkp. (25) 

From the Dyson equation we have 

G{k;z) = i , where Ek{zy.^el + hY,{k; z). (26) 

z - Ek{z) 

Poles of G{k; z) are thus seen to be the isolated (see § A.l) solutions of 



Here we are referring to metallic systems. For non-metallic ground states, this problem of divergent self-energy contributions 
does not arise. Since non-metallic systems are by definition non-Fermi liquids, they do not concern us here. 
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Ek{z) = z, (27) 

which is equivalent to Eq. ( [l^ ) specialised to the uniform isotropic system under consideration. According to our 
conventions, kp has to be identified with kp and we need only to have kp for indicating a k infinitesimally larger 
than kp. As we have mentioned in § IV. C, indeed Ekp{sp) = ep must be satisfied. However, since Yi{kp]z) has a 
branch point sX z ~ ep, it follows that G{kp; z) cannot have a pole ai z — ep, since a branch point is not an isolated 



singularity (see § A.l). Moreover, as we have discussed in § III, Eq. (|14[), or its equivalent, Eq. (27), cannot have a 
complex- valued solution, so that there cannot be any question of a 'pole' of G{k; z) changing the sign of its imaginary 
part. The question, therefore, arises as to how, despite the fact that G{kp;z), or G{kp;e) in Migdal's notation, has 
no pole at z = ep, the above-introduced Migdal (1957) theorem could hold. 

To answer the above question, we note that since, for k = kp and z ep, z — Ek{z) approaches zero (by continuity 
and the fact that Eq. (foh applies) , in principle it may be possible to write 



G '^{kp;z) = ^{z ~ [el^ +nT,{kp;z)]} a{z ~ ep) + L{z), as z ep, (28) 

with L{z) satisfying liuiz^ep L{z)/{z — ep) = 0, that is L{z) = o{z — ep) for z ep. In Eq. (p8|), a stands for a 
constant to be specified below. Equation ( p8| ) implies G~^{kp; z), and thus T,{kp; z), to be at least once continuously 
difFerentiable with respect to z at z = ep (see footnote ^); mth continuous differentiability of G^^{kp; z), and thus 
of Ti{kp; z), with respect to z at z — ep amounts to the condition G~^{kp; z) SjLo '^ji^ ~ + E{z) with {aj} 
finite constants and limz-,eF E{z)/{z — e^)™ = 0. ^ Consequently, dY^{kp;e) / de\e=ep is well-defined and one has 
a = ?i^^(l — hdT,{kp; e)/de\e=£f,) = l/{hZkp) (see Eq. ( p3| ) above). Using the representation z — ep = gex.p{iip), 
with g a non-vanishing constant, to be let to approach zero, from Eq. ( p^ it can straightforwardly be shown that 
Eq. ( p5| ) remains intact, even though z = ep is not a pole (i.e. an isolated singularity) of G{k; z). 

Through the Kramers-Kronig relation (see Appendices B and C) for ReS(fc;e) in terms of ImE(fc;£) and the 
asymptotic relation ImS(fc;e) ^ ^ak{e — ep)"^, for e^ep, it can readily be deduced that (see Appendix C) 

ReT.{k;e)r^j:{k;ep)+Pk{e-ep), with (3k < tor k ^ kp. (29) 



Concerning the sign assigned to for k close to kp, as is seen from Eq. (31) below, a positive (3kp would imply 
Zkp > 1 which is impossible on account of a combination of the following three reasons. 

(i) By definition (see Eq. ( p4[ ) above) < n(fc) < 1; this can also be traced back to the Pauli exclusion principle. 

(ii) According to Eq. (|25|), Zkp = "(/c^) — n(fc^) so that in view of (i), \Zkp\ < 1. 

(iii) Zkp < is excluded owing to ImS(fc;e)^0, for e^^. 

In Appendix C (see text following Eq. (|C9|)) we demonstrate that f3k < 0, for fc in a neighbourhood of kp, directly 
follows from a general sum rule involving ImS(fc; e). We therefore have ^ 



'^^ Because of G{z*) — G'{z), it follows that L{z*) — L*{z). Further, unless e £ (/ijv, /ijv+i), in general lim^|o{L(£ -f 177) — 
L{e — ir])} ^ 0. This is a consequence of the fact that z = /ijv = sf and 2 — ^n+i are branch points of G{k; 2), as discussed in 
§ HI- 

We note that, when a function, say f{z), is analytic inside and on the boundary of a simply-connected region TZ of the 
complex 2-plane, it is infinitely-many times differentiable in TZ, so that once-differentiability of f{z) at a point, say zq, interior 
to TZ implies infinitely-many times differentiability of f{z) at 2 = zq. In the case under consideration, since 2 = ef is a 
branch point, it cannot be interior to any simply-connected region of analyticity of G{kF;z) (see § A.l). Consequently, mth 
differentiability of this function at 2 = ef does not imply its (m -I- l)th differentiability at this point. Consider for instance 
f{z):={z — zq)^ with p > a real non-integer. This function has two branch points, 2 = 20 and I/2 = 0. With m an integer 
satisfying m < p < m + 1, f{z) is seen to be m times, but not m-\- 1 times, differentiable at 2 = zq. 

Consider the analytic function j(z) = u{x, y) 4- iy{x, y), where x -\- iy = z and u{x, y) and v{x, y) axe real. Upon setting y 
equal to zero, we obtain f{x) = u{x, Q)-\-iy{x, 0). Replacing x in the argument of / by 2, we obtain f{z) = u(2, 0) -|-iv(2, 0) (in a 
similar manner one obtains f[z) = u(0, —iz) -|-iv(0, —iz)). This expression prescribes how to obtain E(A;; 2) from the knowledge 
of E(fc;e). In the present case, as 2 ep, we have: E(fc; 2) ~ E(fc;£F) -|- I3k{z — ep) — iak{z — ef)"^, for < arg(2) < tt, and 
E(fc; 2) ~ E(fc; ef) + Pk{z — ef) + iak{z — ef)'^, for — tt < arg(2) < 0. The self-energy as specified in this way is seen to preserve 
the reflection property E(fc; 2*) — E*(fc; 2) for Im(2) 7^ (see Eq. (p^). 
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^{k;s)^^{k;sF)+f3k{e-eF)Tiak{e-eFf, e^ep- (30) 
A simple calculation reveals that for Z^p in Eq. ( p5| ) (or Eq. (|28 

Zkp = {I - hhp)-^ (31) 

holds. A similar analysis based on the Kramers-Kronig relation for the self-energy (see Appendix C) yields that when 
ImE(fc;e)-'| £ — £f|'^ for e — > Ei?, the following hold. 

(i) ReS(fc; e) - T.{k;eF)--'{e -Ef) when 1 < cr < 2. 

(ii) ReE(fe; e) — Yj{k] e f)^' {e — £f) In |e — £_f| when a = 1 ("marginal Fermi liquid" (Varma, et al. 1989, Littlewood 
and Varma 1991, Kothar, et al. 1991) (see Appendix C). 

(iii) ReS(fc;e) — S(fc;ei?)~'|e — when < cr < 1 (see § IV. C, text following Eq. (p3|), and identify ct with 70; 



see also Appendix D and identify a with 1 — 270 in Eqs. (D21) and (D22) and with 1 — 70 in Eqs. (D28) and (D29)) 



It can straightforwardly be shown (using the same strategy as in the case of ct = 2 employed above) that, for 
< (7 < 1, Zkp on the right-hand side of Eq. ( ^5|) is vanishing, whereas for 1 < cr < 2 it takes a finite value. Vanishing 
of Zkp corresponds to disappearance of the Landau QPs on the Fermi surface. Although for 1 < cr < 2, Zkp is 
non- vanishing, it is evident that the smaller the cr, the shorter are the life-times of the QPs close to the Fermi surface, 
if QPs can at all be meaningfully defined (recall condition (B) introduced in § I). We should emphasis that as G{kF\ z) 
does not possess an isolated pole at z = Ef = fJ'N, but a branch point, it follows that even in the case of the Landau 
Fermi liquids, the QPs on the Fermi surface are not truly infinitely long-lived but only so in an asymptotic sense. 

We conclude that what Migdal (1957) describes as 'leaving of a pole from the interior of contour C upon changing k 
from kp to kp' is to be understood as follows: for k = kp, Eq. ( [l^ ) (or Eq. (p7|)) cannot be satisfied when z — ef = fJ^N, 
that is Ei^+{fipf) ^ fipf; rather we have (/xjv+i) = /ijv+i and, as can be viewed from Fig. 2, z = ^n+i indeed does 
not lie in the interior of C. 

In closing this Section, we indicate the following observation. In §§ I and IV. C we argued that Fermi liquids are 
distinguished by two specific aspects, namely that of condition (A), that is continuous differentiability of S(fci?;£) 
with respect to £ in a neighbourhood of £ = ef, and of condition (B), that is continuous differentiability of S(fc; ef) 
with respect to A; in a neighbourhood of k — kp- It is owing to these that the Landau quasi-particle dispersion 
as presented in Eq. ( pO| ) can be rigorously obtained from the quasi-particle equation, Eq. (|l8|), which is applicable 
to all isotropic metallic systems of spin-less fermions; conditions (A) and (B) are prerequisite to assigning a unique 
mass as well as velocity to an elementary excitation at and in a close vicinity of the Fermi surface (see Eq. (|2l| ) 
above). However, as we have observed in the present Section, in obtaining Eq. (|5|), only Eq. (|8|) combined with 
the condition limz^ep L{z) = (which together embody aspect (A)) played a role. Therefore, a non-vanishing Zkp 
which is universally considered as the hallmark for Fermi liquids, is only a necessary condition for a metallic system 
to qualify as a Fermi liquid, but by no means a sufficient condition (see § I). In other words, metallic systems with 
Zkp 7^ can in principle be non-Fermi liquids in the above sense. 

VI. (NON-)FERMI LIQUIDS AND NON-DIFFERENTIABILITY OF THE SELF-ENERGY WITH 
RESPECT TO MOMENTUM CLOSE TO AND ON THE FERMI 'SURFACE' 



(32) 



From Eq. (|24| ) , making use of Eq. ( |26D , one obtains 

dn{k) f dz 1 dEk{z) 



dk Jc 2ni (^z-Ekiz)) dk 



In view of the singular behaviour of n(fc) at fc = fc^?, in employing Eq. (32) we approach kp from either left or right, 
i.e. by 9n(A:)/9fc|j,^j,T we imply left and right derivatives, respectively. 

For our further discussions it is convenient to introduce some notational conventions. We denote the circular contour 
circumscribing £i? in Fig. 2 by Cs (subscript s refers to 'singular', in view of the singularity of G{kF', z) at z — ef) and 
its complement, C\Cs, by C,- (subscript r refers to 'regular'). Consequently, we write n{k) — ns{k) + nr{k) with n^^k) 



Recall that indicates that the corresponding asymptotic relation is correct up to a multiplicative constant. 
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and nr{k) corresponding to the z- integrals along Cg and Cr, respectively (here we assume the radius g characterising 
Cs to be infinitesimally small). 

In arriving at the result in Eq. (p5|), we have made use of nr{kp) — r\r{kp) so that Eq. ( ^ ) in fact amounts to 
ns{kp) — ris{kp) = Z^p- A similar property, namely dnr{k)/dk\i^^^.- = dnr{k)/dk\i^^^.+ , may not be apphcable in 
general, in particular since dnr(k) j dk\^^^^ may be unbounded. Below we demonstrate that: q i) for Fermi liquids 
dn{k) / dk\f^^j,- is either finite, in which case it follows that S(fc; z) / dk\i,^^,- is finite for all z, or infinite, in which case 
dlik/ dk\f,^j,- is also infinite and consequently 9{S(fc; z) — £f)} / dk\f^^j,- is infinite for all z (note that, for Fermi 
liquids, S(A:; sp) is by definition a continuously differentiable function of fc in a neighbourhood of kp)', ii) for non- Fermi 
liquids characterised by Zkp = 0, dns{k)/dk\^,^i^- = 0. Further, although for these systems r\{kp) — n{kp) = holds, 
dn(k)/dk may or may not diverge as k ] kp {k I kp); in the one-dimensional Luttinger model (Luttinger 1963, Mattis 
and Lieb 1965), which is a prototype of one-dimensional systems of interacting spin-less fermions (Haldane 1980, 
1981), with 'weak' interaction — as characterised by a small anomalous dimension, ^ namely < a < 1 — , dn{k)/dk 
does diverge for k —>■ k^, and consequently i9I](fc; z)/9/c|j,^j,- is infinite for all z; here, contrary to the case of Fermi 
liquids, dTilk; ep)/dk\j^^f,- is not necessarily finite. This aspect can be explicitly verified from the available explicit 
expression for the spectral function corresponding to the retarded single-particle Green function (see Appendix D). 

Now we proceed with demonstrating the above statements. In cases where Eq. ( p8[ ) and lim^^j^^ L{z) = are 
satisfied, making use of the Cauchy residue theorem, one obtains 



dn^ik) 



dk 



d'Y,{k-e) 



dns(k) 



c='F ' dk 



k=k 



+ 



= 0. (33) 



In view of the derivative with respect to e, S(fc;e) in Eq. (|3^) may be replaced by I](fc;£) — T,{k;eF). From this 
it follows that a possible divergence of dr\s{k)/dk for k ] kp cannot be due to Yj{k\£p) and therefore a divergent 
dns{k)/dk for k ] kp is not in conflict with the fundamental assumption of the Fermi-liquid theory, namely that of 
continuous differentiability of Y,{k;ep) in a neighbourhood of k = kp. We note in passing that, through employing 
Eqs. (|o|) and (|3l]), we can also write 5ns(fc)/9fc|j,^j,- — dZk/dk\k=kF, where Zk:—{1 — hPk)^^. ^ In Appendix C we 
demonstrate that divergence of dPk/dk (which in view of the latter result implies that of dr\s{k) /dk) at any k = ko, 
e.g. k = kp, signals divergence of dT,{k; z) / dk\k=ko for all z, Im(z) ^ 0; in Appendix D we explicitly show this to 
be the case for the (one-dimensional) Luttinger model for spin-less fermions. With reference to Eq. (^2|), as well as 
Eq. ( |2^ ) where Ek (z) is defined in terms of E(fc; z) , divergence of ^n^ (fc) /dk a.t k = kp is thus seen to imply divergence 
of dn{k)/dk a.t k = kp. This completes demonstration of our point (i) above. It is relevant to mention that for a 
three-dimensional isotropic system of fermions interacting through a hard-core potential, qualifying as a Fermi liquid, 
dn{k)/dk has been shown to be logarithmically divergent for k ^ kp (Belyakov 1961, Sartor and Mahaux 1980). 
^ On the other hand, for a similar system of electrons, interacting through the long-range Coulomb interaction, 
no divergence is observed in dn(k)/dk, k k^, within the framework of the random-phase approximation (RPA) 
(Daniel and Vosko 1960); calculations of n(fc) beyond the RPA are also available (Geldart, Houghton and Vosko 1964), 
however the results do not provide information with regard to regularity or otherwise of dn{k)/dk at k = k^. 

For non-Fermi liquids, such as the marginal Fermi liquid considered in Appendix C, the result 9ns(fc)/9fc|j.^j,^ = 

follows exactly for the same reason that ns{kp) — for these systems: through parametrising Cg in terms of the 
circular-polar angle if, one observes that despite singularity of the integrand at ep, the (^-integration over [— 7r,7r) 
yields an identically- vanishing contribution. As our considerations in § V have made explicit, at the root of this 



Similar considerations can be given for the case k = kj,. 
'"^ The non-interacting case corresponds to a = 270 — 0. See § IV. C and Appendix D. 

We should emphasise that here Zk is merely a generalisation of Zkp and does not necessarily have the same physical 
significance as Zkp does. To appreciate the relevance of this remark, one should realise that for k 7^ kp, Eq. ( [l8| ) may not have 
a real- valued solution. 

In our judgement, nothing precludes the possibility that this divergence is an artifact of the second-order perturbation 
theory employed in the treatment by Belyakov (1961). 

For a calculation of n(A;) at large k as well as some general discussions concerning n(fc) see also Yasuhara and Kawazoe 
(1976). 
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property lies the non-difFerentiability with respect to e of Tj{kp\£) at ep- We emphasise that this is sufEcient for 
disqualifying a system to be a Fermi liquid; our discussions in § IV. C have made evident that this non-differentiability 
condition stands in the way of obtaining the quasi-particle energy dispersion in Eq. ( pO| ) from Eq. (18). 

For a general non-Fernii-liquid metallic system, in arbitrary c?-dimensional spatial space (if indeed such a system 
exists in d > 1), our knowledge with regard to behaviour of dn{k)/dk in some neighbourhoods of k — k"^ cannot 
surpass that of dn[k)/dk corresponding to Fermi-liquid systems. ^ Since, by definition, to a non-Fcrmi-liquid 
metallic system no quasi-particle energy dispersion similar to that in Eq. ( |2C| ) can correspond (see above), and this 
can be effected through a non-continuous differentiable S(fci?;e) in a neighbourhood of £ = ep, the non-continuous- 
differentiability with respect to k of I](fc;ei?) pertaining to these systems cannot in general be ruled in but also not 
ruled out. In d = 1 dimension, the exactly-solvable Luttinger (1963) model provides us with the opportunity to 
arrive at a rigorous statement, however. In this model, for k close to kp one has (see, e.g., Voit 1993b) n(fe) ~ 
1/2 — CiSgn(fc — kp)\k ~ kp]" + C2{k — kp) where Ci and C2 are constants and a = 2^q, the 'anomalous dimension', 
is not universal but depends on the nature and strength of the particle-particle interaction in the system under 
investigation; for the one-dimensional Hubbard model corresponding to a finite repulsive on-site interaction, [/, and 
away from the half-filling of the Hubbard band, < a < 1/8, with the upper limit achieved for infinitely large U. 
One observes that for < a < 1, which in view of the latter observation should amount to a wide range of values for 
a, n(fc), although continuous, is not continuously differentiable in a neighbourhood of fc = kp; dn(k)/dk diverges as 
k ^ k'^ . 



VII. SUMMARY AND CONCLUDING REMARKS 



In this work we have presented a critical analysis of a celebrated theorem due to Luttinger (1961) concerning 
energies and life-times of low-energy QP excitations in interacting systems. This theorem has played the crucial role 
of classifying all metallic systems in spatial dimensions larger than unity as Landau Fermi liquids, with the explicit 
assumption with regard to applicability of the many-body perturbation theory to all orders to these systems. |^ In 
this work we have demonstrated that Luttinger's (1961) proof involves an implicit assumption with regard to the 
dispersion of the QP energies, which assumption we have explicitly shown to be specific to Fermi-liquid systems. We 
have therefore shown that Luttinger's proof amounts to a demonstration of consistency of the mentioned implicit 
assumption with the property Iml](fc; £)~'(e — £f)^, for e ep, to all orders of the many-body perturbation theory. 
It follows that, contrary to the commonly-held view, the many-body perturbation theory does not necessarily break 
down when applied to systems whose self-energies are of non-Fermi-liquid type. In the absence of any theorem to 
replace the Luttinger theorem, one could reasonably conjecture that it is most likely that, in particular in two spatial 
dimensions, non- Fermi-liquid- like metallic systems if not abundant, should not be rare and that these may be correctly 
addressed within the framework of the many-body perturbation theory. In this connection it is important to point 
out the following two observations. 

First, any static and local approximation to the self-energy operator implies a Fermi-liquid behaviour for the 
QPs, no matter how ingenious such approximation may be. The entire body of the energy-band methods based 
on the conventional density-functional theory, involving a local and energy-independent effective potential (for a 
comprehensive review see Dreizler and Gross 1990), pertains to this category of approximation frameworks. The (self- 
consistent) Hartree-Fock scheme which involves the non-local static exchange self-energy, may describe non-Fermi 
liquid metals, however only those whose Z^.p = 1. Second, in evaluating the self-energy operator for a Coulomb 
system, within the framework of the many-body perturbation theory, it has to be realised that in metallic systems the 



Here we employ the notion of 'neighbourhood' in a loose sense, since, as we have mentioned earlier, here by dn{k)/dk, 
k — > kp, we express left or right derivative of n(fc), thus excluding k = kp from the interior of the interval over which n(fc) is 
being differentiated. 

This on account of the fact that the proof of this theorem has been based on an infinite-order perturbation expansion of 
the self-energy operator. 

For uniform systems of electrons interacting through the Coulomb potential, the derivative with respect to k of the Hartree- 
Fock self- energy is logarithmically divergent at k = kp (Ashcroft and Mermin 1981, p. 334) — see Appendix B. Since 
is independent of e (Appendix B), it is therefore analytic over the entire z-plane and thus the associated Zkp = 1 (see § V). 
Therefore a metal within the Hartree-Fock scheme may be non-Fermi liquid solely on account of violation of condition (B) 
introduced in § I. 
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long range of the bare Coulomb interaetion gives rise to divergent self-energy contributions (from the second order 
in the interaction onwards). Thus, for such systems, perturbation expansion has to be in terms of the dynamically- 
screened electron-electron interaction function. In calculating this function, at least all polarisation diagrams of the 
random- phase approximation have to be taken into account, for all these diagrams have unbounded contributions 
and only in combination give rise to a finite result. Our considerations with regard to the problem of unbounded 
self-energy diagrams have led us to draw the additional conclusion that, irrespective of our above-indicated finding, 
Luttinger's (1961) proof cannot have bearing on metallic systems of electrons interacting through the long-range 
Coulomb interaction function. This follows from the fact that Luttinger's proof has been based on a perturbation 
expansion of the self-energy operator in terms of the bare electron-electron interaction. In this connection we point 
out that analytic properties of a series of functions of complex variable do not coincide with those of the constituent 
terms, unless the series be uniformly convergent (Whittaker and Watson 1927, pp. 91 and 92, Titchmarsh 1939, pp. 
95-98). Further, the sum of a series which is not absolutely convergent depends on the order of summation of terms 
in the series (Whittaker and Watson 1927, pp. 18 and 25). It is evident that, in particular, perturbation series which 
involve singular terms cannot be absolutely or uniformly convergent. 

We have further elaborated upon the analytic properties of the single-particle Green function and the self-energy 
operator as functions of a complex-valued energy parameter, z. In particular we have pointed out that the commonly- 
used equation for energies of the QPs in an interacting system can yield, if any, only real-valued solutions. To 
obtain complex- valued solutions (concerning systems in the thermodynamic limit), the self-energy operator has to be 
analytically continued into a non-physical RS of the complex energy plane. 

We have paid careful attention to the precise nature of the singular points of the single-particle Green function and 
the self-energy operator. For instance we have pointed out that the Fermi energy is not an isolated singularity of the 
Green function and therefore cannot be a pole. In the light of this, we have considered a celebrated theorem due to 
Migdal (1957) and discussed how its proof is independent of the commonly-made assumption, that on the Fermi surface 
the Fermi energy were a pole of the single-particle Green function. We have explicitly shown that for the momentum 
distribution function n(fc) pertaining to a metallic system to be discontinuous at k = kp by sl finite amount — the 
magnitude of this discontinuity being, according to the Migdal theorem, equal to Zkp — , it is only necessary that 
the self-energy T,{kF;£) be a continuously differentiable function of £ in a neighbourhood of £ = ep- Consequently 
metallic systems with Zkp ^ are not necessarily Fermi liquids. Fermi liquids, we have shown, have in addition 
the property that their corresponding S](fc;£i?) is a continuously differentiable function of fc in a neighbourhood of 
k = kp- In spite of this, T,{k; e) pertaining to Fermi liquids may not be continuously differentiable function of fc in a 
neighbourhood of k = kp when e ^ ep. These observations, which we have explicitly examined on the self-energies 
of Fermi-, marginal Fermi- and Luttinger-liquids, clearly demonstrate how uncritical Taylor expansions of various 
functions of k and e associated with the self-energy can lead to incorrect description of the physical properties of the 
systems under consideration at low energies. 



ACKNOWLEDGEMENTS 

With appreciation I acknowledge support by the Max-Planck-Gesellschaft, Germany. It is a pleasure for me to thank 
Professor Peter B. Littlewood and members of the Theory of Condensed Matter Group for their kind hospitality at 
Cavendish Laboratory where this work was completed. I record my indebtedness to Girton College, Cambridge, for 
support in the course of my present association with Cavendish Laboratory. I take this opportunity and extend my 
heartfelt appreciations to Dr Martin W. Ennis for generously and patiently sharing his insight especially into the 
music of J.S. Bach and kind hospitality. 

APPENDIX A: MATHEMATICAL PRELIMINARIES 

Since in the present work we repeatedly encounter a number of specific mathematical notions, we devote this 
Appendix to a brief exposition of these. 
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1. On the types of singularity 



A point at which a function g{z) of complex variable z is not analytic ^ is called a singular point (for details of 
what follows see, for example, Whittaker and Watson 1927, pp. 102 and 104, Titchmarsh 1939, pp. 89-95). Such a 
point is either isolated or non-isolated; zq is isolated if there exists a i5 > such that inside the region |z — zo| < 6 
there exists no other singular point of g(z) than zq. Otherwise zo is not isolated. 

A singularity may be removable; such singularity corresponds to a point zq at which g(z) is not defined, but 
lmiz—,zo g{z) exists. Thus z = is a removable singularity of g{z):=sm{z) / z. 

Limiting (or accumulation) point of an infinite sequence of poles, is not classified as a pole and thus is considered 
as an essential singularity, and here, a non-isolated essential singularity. For instance, the sequence of poles of 
g{z)-= J2'^=o ^/('^•[l + a^"z^]), with a > 1, have z = as their limiting point which is not isolated. This function has 
no Taylor or Laurent series expansion (see § A. 2) over any domain of the z-plane which has z = as its interior. 

Let g{z) be a single-valued function throughout a domain "D at whose interior point zq, g(z) is singular. Suppose 
that the principal part of the Laurent series expansion (see § A. 2) of g{z) around z — zq terminates with the term 
a-n/{z — zo)"j with a_„ a non-vanishing constant. In such case, zq is called a pole of order n. Poles are thus by 
definition isolated singularities. 

If the principal part of the Laurent expansion of g{z) around z = zq does not terminate (i.e., if there exists no Hq 
such that a-n = for all n > no), Zq is an isolated essential singularity of g{z). The function g{z):=exp{z) has one 
such singularity at the point of infinity, i.e. at z = 1/C where C = 0- 

Branch points belong to the class of singular points and concern multi-valued functions. Let g{z) be one such 
function. By traversing a closed contour which circumscribes only one branch point of g{z), one obtains, upon 
arriving at the starting point zi, a different value for g{zi), indicating change of the initial branch of g{z) into a 
different branch; for a branch point of order p, the original branch is recovered after completion of p revolutions along 
the mentioned contour. Thus (z — zq)^^^ has a third-order branch point at z = zq. Functions can also possess branch 
points of infinite order; for g(z):=ln(z), z = and 1/z = are such points. 

2. Taylor, Laurent and asymptotic series 

In the Abstract of this work we have referred to the expression ImE(fc;£) ^ =pafc(e — sp)^, e^ep, for e — > ep, as 
an asymptotic relation. Here we specify what asymptotic relations and asymptotic series are and in what essential 
respects these series differ from the Taylor and the Laurent series (Whittaker and Watson 1927, pp. 93, 94 and 100). 
We also indicate the interrelation between branch points (see § A.l) and divergent asymptotic series. 

When g[z) is analytic at z — zq, then, by definition, there exists an open domain 2? of which zq is an interior 
and over which g(z) is analytic. Within a circle around zq embedded within T), g{z) can be represented in terms 
of a Taylor series: g{z) — X]tiLo'^"('^ ^ ^o)", with unique coefficients a„. Thus {(z — zo)"|n — 0,1,...} can be 
considered as a complete basis in a (finite) neighbourhood of zq for any function that is analytic at z — zq (note 
that this set does not contain such term as, for instance, (z — zg)'^ with a non-integer). This basis can be extended 
to {(z — Zo)"|n — — m, — m + 1, . . . , 0, 1, . . .} in order to form a complete basis in a (finite) neighbourhood of zq for 
representing any function which has a pole of order m or lower at zq and is analytic in a neighbourhood of z = zq; 
the resulting representation is the well-known Laurent series expansion (§ A.l). It cannot be extended, through 
increasing m in {(z — zo)"|n = —to, — m + 1, . . . , 0, 1, . . .}, for representing functions that possess a branch point at 
z — Zq (see § A. 3). It is important to point out that the Taylor and the Laurent series expansions are very special 
in that they provide uniform representations of an analytic function; when these series are convergent for any g and 
ip in z — Zq = gexp{iip), they are uniformly convergent for all values of ip, ip G [0, 27r). Branch points are peculiar in 
that they do not allow for any uniform representation (note the '=f' in the above expression for ImS(fc;e)). We shall 
return to this point further on in this Section. 

Asymptotic series (in the sense of Poincare) (Whittaker and Watson 1927, Ch. VIII, Copson 1965, Dingle 1973, 
Lauwerier 1977) are more general than the Taylor and the Laurent series in that they are in terms of "basis functions" 
that are not necessarily of the form (z — Zq)", with n possibly negative, zero and positive integer. To an asymptotic 
series corresponds a so-called asymptotic sequence: The set {(^„(z)|n = 0,1,...} is an asymptotic sequence with 



Analytic, regular and holomorphic are alternative but equivalent designations. 
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respect to zq provided it possesses the property 0„+i(z)/(/)„(z) — > 0, denoted by 0„+i(z) = o((^„(z)), for z zq. 
Thus {(z — zo)"|'^ = —m, —m+ 1, . . . , 0, 1, . . .} is an asymptotic sequence with respect to zg. We point out that some 
authors, for example Whittaker and Watson (1927, Ch. VIII), reserve the designation 'asymptotic series' for those 
series that are both asymptotic (in the above sense) and divergent. We do not foUow this restrictive convention. 

As we have mentioned above, an asymptotic series expansion of a function may not be convergent. It may not 
also be uniform; functions with non-uniform asymptotic series expansions for z ~> zq are characterised by possessing 
different asymptotic series for different sectors of the z-plane around z = zq. Thus, for instance, 17(2):= exp(z) + 
exp(— z) tanh(l/z) has the following asymptotic forms for z ^ 0: g{z) ~ 2cosh(z) ~ 2 + z^ + . . ., when Re(z) > 0, 
and g{z) ~ 2sinh(z) ^ 2z + z^/3 + . . ., when Re(z) < (Lauwerier 1977, p. 11). The fact that in different sectors 
around a point in the complex z-plane a function can have different asymptotic expansions with respect to the same 
asymptotic sequence, is referred to as the Stokes phenomenon (Watson 1952, Dingle 1973); 'Stokes lines' are thus the 
branch cuts in our discussions. 

The Taylor series of analytic functions based on point zq are convergent uniform asymptotic series for z — > zq 
(similarly for the Laurent series). It can be shown (Lauwerier 1977, pp. 12-14) that an analytic function can be 
associated with even a divergent asymptotic series corresponding to the asymptotic sequence {(z — zo)"|n 0, 1, . . .}; 
the given (divergent) series is its asymptotic expansion. Thus, for instance, after Borel (or Euler) transformation 
(Whittaker and Watson 1927, pp. 154 and 155, Dingle 1973, pp. 405-408, Lauwerier 1977, pp. 45-50) of the 
divergent series f{z):— X]i^o(~^)""''^"' '^^^ formally obtains /b(z):= X]^o(~1)"-^" = (1 + ^)^^- Through the Borel 
back-transformation of Jb^z), f b{z)-= /o°° exp(--g;)/B(xz), one obtains a function, i.e. f^iz), which is analytic 
in the sector —it < arg(z) < tt of the z-plane. |^ Through replacing fsixz) by its formal geometric expansion 
1 — XZ+ (xz)'^ — • ■ •, the tcrm-by-term evaluation of the latter integral yields the original divergent series. We observe 
that the divergence of the original asymptotic series is closely associated with the restricted sector of the z-plane 
around z = over which f b{z) is analytic. 



3. Many- valued functions: Physical and non-physical Riemann sheets 

Many-valuedness (Whittaker and Watson 1927, pp. 96-98, Titchmarsh 1939, pp. 138-164) is a generic property of 
correlation functions of complex energy pertaining to systems in the thermodynamic limit (see § A. 4). In particular, 
the single-particle Green function and the self-energy operator are many- valued (see § II). 

An n-valued function of complex variable z over domain V may be thought of as embodying n branches of a single- 
valued function over the extended domain consisting of the union of n replicas of P; since these domains signify the 
same region of the complex plane, they are to be distinguished by considering them as being located on different sheets, 
namely the Riemann sfieets (RSs), of the complex plane. In the present paper we single out one particular branch 
of the many-valued functions that we encounter, such as the single-particle Green function, because of its physical 
significance. To this particular branch we refer as the function in question on the physical RS or the 'physical' function 
(compare with G(e) in Eq. (|j) which we have designated as the 'physical' single-particle Green function). 

Throughout this paper we employ the following notational conventions. Let /(e) be a function of real energy variable 
£, which we consider to describe some dynamical physical property of the many-particle system under consideration. 
By /(z) we denote the analytic function that uniformly approaches /(e) when z ^ e; the so-called 'edge-of-the- wedge' 

theorem (Streater and Wightman 1964) asserts uniqueness of /(z). Q We refer to /(z) as /(z) on the physical RS; 

/(z) denotes the above-mentioned single-valued function over the union of all the RSs. We recall that in the main 
text we have referred to, for example, G(z) in Eq. (^) as the single-particle Green function on the physical RS. 

An example will clarify this. Let /(e):=ln(e), for e > 0. We can identify ln(z) by /(z) since for z — > e, with e > 0, 
indeed ln(z) — > ln(e). One can easily verify that Ln„(z):= In |z| + i{arg(z) + 27rn}, with n an integer (positive, zero 



It can be shown that fsi^) = exp(l/2)i5i(l/z)/z where Ei{z) stands for the exponential-integral function (see, e.g., 
Abramowitz and Stegun 1972, p. 228). We note that exp(l/2:) has an isolated essential singularity at z = (see § A.l) and 
that Ei{l/z) has branch points at 2 = and 1/z = 0. 

^® The condition of uniformity can be shown, a posteriori, to be redundant: the existence of the limit f{z) fi^), as z e, 
with £ varying over an open interval where /(e) is continuous, implies uniformity of f{z) /(e); see p. 75 in Streater and 
Wightman (1964). 
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and negative) and — tt < arg(z) < tt, is the analytic function (i.e. f{z)) of which /(z):=hi(z) is the physical branch; 

thus this branch corresponds to n = 0. The physical branch f{z), like any other branch of f{z), is a many-valued 
function (branch points are not removable [see § A.l] singularities). It has two branch points, one at z = and the 
other at 1/z = 0. Our above convention — tt < arg(z) < tt specifies the branch cut of f{z), which connects the two 

branch points of f{z), to be along the negative e-axis. One can determine all branches of f{z) from the knowledge 

of any of its branches. Two of infinitely many branches of f{z), fi{z) and f2{z) say, that can be obtained by means 
of a direct analytic continuation of f{z), through its branch cut, into two diff'erent non-physical RSs are uniquely 
determined from the following requirements: lim^|o{/i(£ + ^??) — /(£ — ir])} — and lim^|o{/2(£ — *??) — /(e + *??)} — 0, 
for £ < 0. It is not difficult to verify that fi{z) = Ln_i(z) and f2{z) = Ln+i(2;). 



4. A physically-motivated example 

Here we apply the above concepts to a simple model function which accommodates a number of salient features 
of the physical functions which we deal with in the main part of this paper. We shall particularly emphasise the 
role played by the process of evaluating the thermodynamic limit in modifying the nature of the singular points of 
dynamic correlation functions. 

Consider /(e; N, f2):=^ ^ ^o) ^ ©(££ — ~ s); with e^+i > e^. Here cq and ei are finite constants for 

which we assume cq < eg < ei for some values of ^; and f2 indicate that / is a function of the number N of particles 
as well as the volume of the system. In the 'thermodynamic limit' {N — > oo, f2 — > c» and finite concentration 
C:=N/^l), f is only a function of C. A brief glance at the contents of § II should clarify our present choice for /. 

Suppose that in the thermodynamic limit {si+m ~ £i) ~* for any finite value of M . This condition implies that in 
taking the limit, the function /(e; TV, VL) becomes ill-defined for any real e in the interval [eo, ei]; an e in this interval 
will be "pinched" (Itzykson and Zuber 1988, pp. 302 and 303) by poles of f{e]N,Vl). To avoid this problem, the 
thermodynamic limit has to be effected after replacing the real energy variable e by a complex energy variable z, since 
a complex z cannot be "pinched" by the real poles of /(z; A^, Vt) as A^ and Q. approach infinity. 

To be specific, let us now assume that in the thermodynamic limit poles of f{z\ N, 57) populate the interval [eg, ei] 
with a constant density equal to yl x 57. We write 



17 -^-^ Ef — z 

^ A f ' -^ = :J(z;C), for A^ ^ oo, 57 ^ cx). (Al) 

One trivially obtains 

f{z; C) = A{\n{z - ei) - ln(z - eo)}. (A2) 

Evidently, the process of evaluating the thermodynamic limit has led to a dramatic change in f{z; N, 57), transforming 
it into f{z;C) which analytically is distinct from /(z;A^, 57). For instance, f{z;C) has no poles, rather it possesses 
two infinite-order branch points, at z — cq and z = e\. Consequently, contrary to /(z; A^, 57), which is single-valued, 
/(z; C) is a many-valued function of z. 

We shall now explicitly demonstrate that /(z; C) as presented in Eq. ( |A2| ) is the 'physical' branch of /(z) (see § A. 3). 
To this end, consider f{e±iri; C) = A J^^ de' / {e' — e^irj). For 77 | we have l/{e' — e^irj) = V{l/{e' — e)}±inS{£' — e), 
where V denotes the Cauchy 'principal value'. It is then trivially seen that /(e ± 177; C) — A{\n |ei — e| — In |eo — e|} ± 
iTrAQ{e — eo)©(ei — e) for 77 J, 0. Hence, for e < cq and e > ei the 'physical' branch f{e ± irj; C) must be real-valued, 
while for e G (eo, ei) it must satisfy lim,,|o{/(£ + i'7; C) — /(e — irj] C)} — 2T:iA; thus the interval [eo, ei] constitutes 



the branch cut of f{z;C). These conditions are exactly fulfilled by the expression in Eq. (A2); note that ln(z), as 



distinct from Ln„(z) for n ^ 0, stands for the principal-value logarithm (see, e.g., Abramowitz and Stegun 1972, p. 



67). Thus /(z; C) in Eq. (A2) is indeed the analytic continuation of /(e; C) into the physical RS. 

It is instructive to consider two examples concerning continuations of /(z; C) into non-physical RSs. Let 
g(z):=^{Lni(z — ei) — ln(z — eo)}. It is easily verified that for eo < e < ei, lim^|o{/(£ + C) — di^ — iv)} = 0; 
which implies that g{z) is the analytic continuation of /(z; C) from the upper-half plane of the physical RS through 
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the branch cut [eo,ei] into the lower-half plane of a non-physical RS. That is 5(2) is a branch of /(z;C) on a non- 
physical RS. As for the second example, consider h(z):—A{LTi^i{z — ei) — ln(2: — eo)}. For e G (eo,ei) we have 
lim^j^o{/(£ ~ *'7; C') ~ + = so that /i(2) is the analytic continuation of /(z; C) from the lower-half plane 
through the branch cut [eg, ei] into the upper-half plane of a non-physical RS. 



APPENDIX B: ASYMPTOTIC BEHAVIOUR OF THE SELF-ENERGY AT LARGE ENERGIES 



Here we demonstrate that, for 



00, ^(z) ^ S , the Hartree-Fock self-energy. 



I]^^(r, r') = E^(r, r') + S^(r, r') 



-VH{r; [n])S{r 



2h 



v{r — r')g{r, r') 



(Bl) 



where Wff(r; J d'^r' u(r — r')n(r') stands for the Hartree potential, which is a functional of the electronic number 

density n in the ground state, and g for the reduced single-particle density matrix. We note that g in the definition for 
pertains to the fully interacting system and is distinct from the Slater- Fock reduced density matrix go; contrary 
to the former, the latter is idempotent, i.e. gogo = go, whereas gg ^ g. This difference may have some far-reaching 
consequences. For instance, in a uniform isotropic system of electrons interacting through the long-range Coulomb 
interaction function, for S^(fc) evaluated in terms of go, which we denote by (fc), one has (Ashcroft and Mermin 
1981, p. 334) [] 



Sf(fc) 



2e^ 



kpTik/kp), H^):=: 



4.T 



■In 



1 



(B2) 



the first derivative with respect to k of (fc) is seen to be logarithmically divergent at k — kp, an aspect which may 
be an artifact of Sf'(fc) and may not be shared by E-^(fc). 

From the Dyson equation we have T,{e) = Go^ie) — G~^{e). Therefore to deduce the asymptotic expansion of E(e) 
for |e| 00, we need to determine those for Gg "^(e) and G~^{e). The asymptotic series for Go^{e) and G~^{e) can 
be deduced from those pertaining to Go{£) and G{e), respectively, through rehance on the following result from the 
theory of asymptotic analysis (see Copson 1965, pp. 8 and 9): Let f{z) have the following asymptotic expansion 
for \z\ — > 00, where z is a complex variable: f{z) ^ fo + fi/z + -t- . . ., with /o, /i, ... constants, independent 

of z. Then, provided that fo ^ 0, l//(z) ^ I//0 + fi/z + + . . . holds for \z\ — > 00, where /i — —fi/f§, 

h = if? - /o/2)//o^ etc. 

From the Lehmann representation for G(e) in Eq. (|l|) it follows that 



Gis) 



G. 



for 



(B3) 



where 



Goo,(r,r') = A,(r)A:(r') ^ M(r - r'), 

S 

Goo.(r,r') = n^e.A,(r)A:(r') ^ n(s<(r, r') + S>(r, r')) ; 



(B4) 
(B5) 



for S< and 5> see Eqs. ( |B9| ) and (BIO) below. A similar expression to Eq. ( |B3| ) holds for Go(£). Therefore from our 
above considerations with regard to the asymptotic series for l//(z) in terms of the coefficients of that for /(z), it 
follows that S(e) possesses the following asymptotic series expansion 



E(£) 



for \e\ 



(B6) 



^"^ For this system, E^(A;) is divergent but cancels an equally-divergent contribution due to the field of the positively-charged 
uniform background. 

When fo — and /i 7^ 0, then one should apply this result to g{z):—zf{z); from the asymptotic series for l/g{z) that for 
17/(2) is obtained through a multiplication by z. 
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Here we are interested in the leading asymptotic term Soon which our above considerations indicate to be determined 
from Go;oo2 and Gooa as fohows 



For the many-body Hamihonian of the form 



H 



2mp 



u(r) 



7?(r) + - / d'^rd'^r' i!\r)i^^r')v{r - r')^?(r')^(r), 



(B7) 



(B8) 



with 7i(r) the local external potential, making use of the definition for the Lehmann amplitudes and energies as given 
in Eq. (||) and the fact that Em,s\'^ ai.s) = -ff|^'A/.s), one readily obtains 



S<(r,r'):-5]0(M-es)e.A,(r)A:(r') 

S 

= -{^N,o\$\r') [i/,?A(r)]_ |*A,,o), 
S>(r,r'):=5]0(e.-/i)e.A,(r)A:(r') 

S 



(B9) 



(BIO) 



By applying the canonical anti-commutation relations for the field operators in the Schrodinger picture one arrives at 
the following result 



0002(1", I"') = ^ 



2m, 



-V"^ + u{r) +VHir; [n]) 



5{r — r') — -^v{r — r')g{r, r') 



(Bll) 



The corresponding expression for the non-interacting Green function follows from this expression by setting the 
coupling constant of the electron-electron interaction equal to zero, 



G0;oO2(l", I"') = ^ 



2mp 



u(r) 



Sir - r') 



From Eqs. (B7), (Bll), (B12) and (Bl) we obtain the following result 



^ooo = S 



HF 



(B12) 



(B13) 



which completes our demonstration. We note that Eqs. (B3) and (B6) hold equally for G{z) and ^(2;), respectively, 
with merely e replaced by z. 



APPENDIX C: THE KRAMERS-KRONIG RELATIONS FOR THE SELF-ENERGY; APPLICATIONS TO 

ISOTROPIC FERMI AND MARGINAL-FERMI LIQUIDS 

The analyticity of a function of the complex variable z — x + iy in domain V of the z-plane implies that in V 
its real (imaginary) part is up to a constant uniquefy determined by its imaginary (real) part (Morse and Feshbach 
1953, pp. 356-358). For functions, such as f{z), whose singularities are located along the x-axis and whose real and 
imaginary parts, u{x,y) and v{x,y), approach zero not slower than 1/gf, with p > 0, as q:—{x^ + y^)^/^ — > 00, the 
Kramers-Kronig relations establish the mentioned unique relationships between u(x,0^) and v(a;,0^). 

The self-energy Tj{z) is analytic over the entire z-plane with the exception of the real axis. Q Two aspects should 
be addressed before a Kramers-Kronig pair of relations can be set up for the self-energy. First, as the considerations 
in Appendix B have demonstrated, T,{z) ^ Y,^^ for \z\ — > 00, so that we need first to introduce 



This follows from the Dyson equation 'S{z) — G^^ (z) — G ^(z) and the fact that neither Go{z) nor G{z) possesses zero 
eigenvalues for Im(2:) 7^ (see last paragraph in § II). 
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S4z):=S(z)-S 



HF 



(CI) 



which, foUowing Eq. (B6), approaches zero according to ^ Hoo-^/z. Second, as Yi{z) is an operator, it cannot be 
assigned real and imaginary parts; moreover, as 'Si{z) is not gauge invariant, real and imaginary parts of the matrix 
elements {(a|E(z)|/3)}, for a given representation, are dependent upon the choice of the gauge. We need therefore 
generalise for operators the notions 'real part' and 'imaginary part' before the corresponding Kramers-Kronig relations 
can be meaningfully defined. Wc introduce 



S'(z):=l{S(z) + St(z)}, I]"(z):=l{E(z) - Y.\z)}. 



(C2) 



We similarly introduce Ti'^{z) and ^"(z); ^'(z) and zE"(z) {Yj'^{z) and iYj'l{z)) are Hermitian and anti-Hermitian 
components of ^(z) (I]*(z)) and are unique. This uniqueness enables us to deduce from 



S.(£±z?;) = ±—V 



de' 



(C3) 



which follows from the application of the Cauchy theorem together with the above-indicated analytic and asymptotic 
properties of 'E-^{z), the following pair of representation-free Kramers-Kronig relations 



^ J-oo e'-e 

In obtaining these relations we have made use of our convention in Eq. (^ 
S"(£)^S:'(£). 



We note that S'(e) = S 



(C4) 
(C5) 
E^e) and 



1. Fermi liquids 



Now we employ Eq. (C4) and obtain from ^ E"(fc;e) = ImE(fc;e) ~ Tc^fel^ ^ £f)^, £<£F) the associated E'(fc;£); 
in this way we obtain Eq. ( p9|) as well as an explicit expression for [3k in terms of ak and E"(fc;e). To this end, we 
subdivide the interval of the e'-integration in Eq. (C4) into three subintervals, (—00, — A], (—A, -I- A] and (+A, -l-oo), 
where we assume A > \5e\ with 5e:—£ — ep- Upon change of variables and transforming integrals over the negative 
e'-axis into those over the positive e'-axis we obtain 



-2akSe 



r 



de' 



-/2 



fc2 



de 



S"(fc;e' + eF) S"(/c; -c' + ei.) ' 



e' — fc 



(C6) 



The first integral on the right-hand side of Eq. (C6) is standard (Gradshteyn and Ryzhik 1965, p. 59) and one has 
V /o^de' e'^/(e'^ — Se^) = A — ((5e/2)ln|(A -I- fe)/(A — fc)|. Since, in the second integral on the right-hand side 
of Eq. (C6), e' > A and since A > |fe|, we replace l/(e' — 5e) and l/(e' -I- 5e) by their respective geometric series 
expansions, in powers of Se/e', and thus obtain a powers series in Se for the second integral. Collecting terms of the 
zeroth and the first order in Se, while making use of Eq. (29), wc obtain 



E(fc;ej.) = E'(fc;e^^) - J:"'\k) + ^^^(fc), f3k = S\^>ik) 



2Aafc 



(C7) 



where 



1 f°° de' r -1 
Si^Hky.=- _^|s"(fc;e' + e^) + (-l)"S"(fc;-e' + e^)}. 



(C8) 



Note that sf = Mjv, Ofe > and that this asymptotic relation is valid for fc in a neighbourhood of fef, say for < k^Skp 
— see footnote |^ 
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In Eq. (jp7\), signifies the approximate nature of the results in so far as A is finite and, moreover, in evaluating the 
e'-integral over (— A,+A] we have employed a truncated asymptotic expansion for T,"{k;e). As should be evident, 
our choice for a finite A has it root in the requirement A > \Se\. Since Eqs. (|C7|) and (|Cq) are valid independently 



of the value of 6e, it is possible to let A in these equations approach zero; because ±e' + EfWs' for s' -> 0, 



(fc), which involves exists for m = 0, 1, and is equal to Sm'{k)] on the other hand, for m > 2, Sin'{k) 



(0), 



,(A), 



does not have a limit for A | 0. Thus we can write 



(C9) 



We point out that, since, for e' > 0, S"(fc; e' + e^) < and i;"(fc; -e' + £f) > 0, it follows that sin^ (fc) < for all 
odd values of m. Specifically, following Eq. (|C9|), we have f3k < (concerning the range of k over which this result 
can be relied upon, see footnotes ^ and ^0|). As we have indicated in § V (see text following Eq. (p9|)), for fc — > /cf 
this inequality guarantees Zkp < 1. 



It is interesting to note that sif'\k), m 





1 

TT 


Im 






1 

TT 






Im 


■Jfj,+A 



0, 1, can be written in the following appealing form 
dz 



{z - fl) 

dz 
{z ~ fl) 



{E(fc; z) + I](fc; 2fi-z)- 2T,"^{k)} 

2 {nk■,z)~m;2^i~z)} 



(CIO) 
(Cll) 



where the term 21]^^ (fc), which is real- valued, has been introduced in order to render the corresponding integral 
existent (see Eq. (B13)). In Eqs. (CIO) and (Cll) wc have further restored /i which originates from Eq. (C4) and 
which is not exactly equal to ep (see § IV. A). We point out that in these expressions, contrary to that in Eq. (|C^), 
A cannot be set equal to zero, since for A = the integrals enclosed by square brackets (logarithmically) diverge as 
a consequence of the integrands being proportional to l/(z — /i) for 2 — > /i; since, however, the singular contributions 

are real-valued, A can be made as small as desired, provided that it is kept non-zero; thus Sm \k)j m ~ 0,1, is 
well-defined. 

The relevance of Eqs. ( plO| ) and ( |Clll ) rests in the following: The fact that I](fc; z) is analytic over the entire z-plane 
(excluding the real axis outside {iJ.N,fJ'N+i )), implies that if d^S^S\k) / dk'^ , m = 0, 1 (including t = 0), is singular at 
some fc, this singularity must be shared by the integrand of d^Sm \k)/dk^, m — 0,1, for all z, Im(z) 0. In other 
words, Eqs. (CIC) and ( pll ) make explicit that possible singularities in d^Sm \k)/dk^, m = 0, 1, are not specific to 
E(fc; z) for some z — zq, but to S(fc; z) for all z, Im(z) 7^ 0. 

In order to demonstrate the utility of our above considerations, let us assume Y,^^{k), similar to S]^^(fc) (see 
Eq. (B2)), has a divergent first-order derivative with respect to fc at fc = fci? (see Appendix B). For the system 
under consideration to be a Fermi liquid, 'E{k;£p) must, by definition, be a continuously differentiable function of 
k at k = kp (see §§ I and IV. C). From Eq. (C9) it follows that the assumed singularity due to dT,^^{k)/dk at 



®^ Violation of these inequalities implies breakdown of causality and instability of the ground state of the system (see Eq. ( |2^ ) 
in comparison with Eq. (0))^ 

®^ Consider g{k):= dz f{k;z), where the contour of integration C lies inside an open domain V of the complex z-plane 
where f{k; z) is analytic. Suppose that the £th-order derivative (including I = Q) with respect to k of g{k) is divergent at 
k = ko. This divergence cannot be due to some isolated singularity (singularities) of f{k; z) / dk^\k=kQ along C, for analyticity 
of f{k; z) enables one to deform C inside T>, thus avoiding the mentioned singularity (singularities), without changing the value 
of the integral (Cauchy's theorem); exceptions to this concern the cases where d^f{k;z)/dk^\k^ko is diverge nt ei ther a t the 
end-points of C, which cannot be dislocated, or at two neighbouring points which 'pinch' C. Note that in Eqs. (CIO) and ( |Cll| ) 
the contours of integration are arbitrary, as long as they connect ^ + A with the point of infinity (the 'end-points') and are 
located on the upper-half of the z-plane (signifying T>). Therefore, excluding end-point and 'pinch' singularities, divergence 
of g{k) /dk''' at k = ko implies divergence of f{k; z) /dk^ at k = ko for all z inside T>. In Eqs. ( CIO ) and (Cll) one of 
the end-points, namely the point of infinity, is harmless, as the integrands in both expressions are vanishing at this point; by 
choosing ^ = (/ijv + ^jv+i)/2 and A = (/ijv+i — ij.n)/2, owing to the assumption of continuous-differentiability with respect to 

k of S(fc; /ijv) and E(fe; fJ-N+i) in a neighbourhood of fc = kp, a possible divergence of Sm \k) or dSm \k)/dk, m = 0, 1, at 
k — kp cannot be ascribed to an end-point singularity. 
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k ^ kp must therefore be compensated by a counter contribution arising from With reference to our above 

discussions, in such an event 9{S(fc; z) + E(fc; 2fi — z) — 21]^^ (k)} /dk must be singular aX k = kp for all z, lm{z) 
(see specifically footnote ^). For the special choice of z = /i + A, with ^ = {^n + /xjv+i)/2 and A = {^n+i — A'w)/2, 
z) and S(fc; 2/i — z) are continuously differentiablc functions of fc in a neighbourhood of k ^ kp (by the Fermi- 
liquid assumption); nonetheless, in accordance with our expectation 9{S(fc; z) + S(fc; 2/Lt — z) — 2J^^^{k)}/dk remains 
divergent aX k = kp owing to the assumed divergence of dYj^^{k)/dk ai k — kp . On the other hand, for an arbitrary 
z the continuous differentiability with respect to k of S(fc; z) and ^(z; 2/i — z) is not guaranteed, so that divergence 
of a{S(/c; z) + E(fc; 2^i - z) - 2E"^{k)}/dk at k ^ kp can be in part due to that of d{E{k; z) + T.{k] 2^ - z)}/dk 
at k = kp. In view of the asymptotic relation in Eq. ( |30| ) (see footnote |4[), it follows that (3k and at may not be 
necessarily continuously differentiablc at k = kp. The divergence of d(3k/dk as k ^ k^ gives rise to a divergent 
dns{k)/dk, as well as a divergent dn{k)/dk, at = fc J (for a discussion of this case see § VI). Since for Fermi liquids 
(here as characterised by Zkp ^ 0), a divergent dn{k)/dk, at k = kp, has proved possible (Belyakov 1961, Sartor and 
Mahaux 1980), we observe that in general even for Fermi liquids S(fc;e) is not a continuously differentiablc function 
of k (in a neighbourhood of fc = kp) for e =/= ej; at £ = ej, f3k and ak are excluded from contributing to S(fc; e) (see 
Eq. (S)). 

2. Marginal Fermi liquids 

Having considered the case of isotropic Fermi liquids in considerable detail, below we briefly deal with the case of 
the isotropic marginal Fermi liquids. 

For marginal Fermi liquids we have Y,"{k;e) = lmJ^{k;e) ^ ~ak{s — ep), as £ ^ ep, with > (see § V). For 
the S^fc; e) corresponding to this S"(fc; e) one obtains an expression which except for the first term is identical with 
that presented in Eq. ( |C6| ). For this first term, which we denote by S^.j^(fc;e), we have (as in the Fermi-liquid case, 
below A > \Se\) 

K,ik;s):=^vj\e' . { ln(A) - In - ^ .s^}; (C12) 

(see Gradshteyn and Ryzhik 1965, p. 59). One observes that here, contrary to the Fermi-liquid case, the limit A J. 
cannot be taken. This is related to the fact that here S"(fc;±e' -f- ep)'^'e' for e' 0, so that only SQ'^\k) has a 
limit for A | 0. Therefore, in the present case, only the first expression in Eq. (|C9| ) is meaningful. Further, from 



Eq. (C12) we observe that the leading asymptotic contribution to S^(/c; e), as e ^ e^, involves In |fc| which is singular 
for 6e — 0. Combining the above results, we have S(fc; e) ^ T,{k;ep) + (3k{£ — ef) In 1^ — £f| + [ik + ictk] {£ — £f), as 
e ^ ep, where 

Pk ^ —, (C13) 

TT 

7fc = -^ln(A) + 5l^)(fc). (CM) 

TT 

The first term on the right-hand side of Eq. (|C14| ) cancels the logarithmically divergent contribution arising from the 
second term as A | 0. It can be readily verified that to the order in which —ak{s — ep) is an exact representation of 
S"(fc; e) for e — + ep, djk/dA = 0. This invariance property correctly reflects the fact that 7^ must not depend on the 
cut-off energy A. 

APPENDIX D: SINGLE-PARTICLE GREEN FUNCTION AND SELF-ENERGY OF THE LUTTINGER 

MODEL; SOME ASYMPTOTIC EXPRESSIONS 

In this Appendix we consider the single-particle Green function and the self-energy pertaining to the Luttinger 
model (Luttinger 1963, Mattis and Lieb 1965) for spin-less fermions in some detail. From the expression for the 



We note that an explicit expression for the single-particle Green function for the Luttinger model with 5-function interaction 
in real space (this model is equivalent with the 1-1-1 Thirring (1958) model, dealt with by Johnson (1961)) has been evaluated 
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former function we can explicitly test the validity of one of our findings in Appendix C, namely that when 9E(fc; z) / dk 
is divergent at fc = fco and z = zq, it is divergent at k = kg for all z. Further, we consider the asymptotic behaviour 
of both G{k;e) and S(fc;e) for e in the close vicinity of ep. In doing so we separately deal with the case where 
k ^ kp and where k ^ kp- In order to remain close to the available literature concerning the Luttinger model, in 
this Appendix we adopt the commonly-used notations and units and therefore in these deviate from the other parts 
of the present work. We specify the new notations as we proceed through this Appendix. In this Appendix h ~ 1. 



1. The single- particle Green function and its derivative with respect to momentum 

We first construct the single-particle Green function G from the spectral function p corresponding to the retarded 
part of the Green function. For this spectral function wc have (Voit 1993b, 1993a) 

Priq,uj) = ■^—^^—^e{i^ + rvoq)Qii^-rvQq)-f(^jo,-^{uj + vorq)j 
( A \7o-i / A \ 

+ (w^-cj,q^-g), for 70^0, (Dl) 
/9r(<?, w) = Po;r(g,w):=='5(tJ - worg), for 70 = 0. (D2) 

Here r = =p specifies the left and right branches respectively of the single-particle spectrum in the Luttinger model, 
to = e — n and k = q -\- rkp . It is important to note that contrary to the 2D and 3D cases, here k and q take on 
both positive and negative values. Moreover q is measured with respect to ^ rkp; k, on the other hand, is measured 
with respect to the origin. Further, A stands for the (finite) cut-off on the range of the interaction in the momentum 
space — the precise value of A is not of relevance here (Solyom 1979; see §§ 9 and 2 herein); T{z) stands for the 
Gamma function (Abramowitz and Stegun 1972, p. 255), 70 = a/2, with a the 'anomalous dimension' (see § IV. C) 
and 7(0, z) denotes the incomplete Gamma function (Abramowitz and Stegun 1972, p. 260). Unless we explicitly 
indicate otherwise, below < a < 1 and thus < 70 < 1/2. With 

Priq,L^y.= --lmG^iq,Lu), (D3) 

TT 

where the denotes the retarded Green function, making use of the fact that lmGr{q,uj) — sgn(a;) ImG^(g,a;), we 
have 

Griq,co) = r dco' Pf:""'^ j 0. (D4) 

J~oo uj-uj' + iri sgn(w') 

From Eq. (|26|), introducing Lo'^.^j-.^rvoq, we obtain (compare with Eq. (^3[)) 

dGr{q,z) ~2 ( duj^.q dJ^r{q,z) 



= G 



dq ^ ' I dq 

Since Gr{q, z) is an entire function of z for all z, Im(z) 7^ (see § II) and, since duj^.^jdq = rvQ is finite, any possible 
divergence of dGr{q, z) / dq, as function of q and z (Im(z) ^ 0), must be due to a divergence of dY,r{q, z) / dq on the 



by Theumann (1967). The corresponding expression in Eq. (2.19)t (here subscript T denotes an equation from Theumann's 
work) turns out to be identically vanishing for the case of non-interacting Luttinger model. To verify this, note that the coupling 
constant A of interaction as presented in Eq. (1.4)t, determines Q. in Eq. (2.3)t which in turn fixes /3 in Eq. (2.14)t. It is seen 
that for non-interacting particles, corresponding to A = and thus /3 = 1, Gi{k,u)) is identically vanishing, owing to sin7r/3 in 
Eq. (2.19)t. We have not attempted to identify the origin of this shortcoming of G\{k,uj) in Eq. (2.19)t. 
®* Thus a better notation for q would be qr in order to emphasise this r-dependence of the origin. 

®^ It is readily verified that Gr{q, z) = duj' Pr{q, uj')/{z — u}'), for Im{z) 7^ 0; in view of Eq. (|^, Gr{q, io) — lim,,j^o Gr(g, (jJ± 
irj), for Lo^O, which can be easily shown to coincide with that presented in Eq. (D4). We point out that, with po-r{q,(^) = 



S{uj — ojr-q) (see Eq. ( |D2| )), one obtains Go;r{q,i^) = ~ i^riq + sgn(cij,'!.g)] and consequently Go;r{q,z) = l/[z — cur-q], for 
Im(2) j^'o. 
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right-hand side of Eq. (p5|). In order to be able to evaluate dGr{q, z)/dq, which according to Eq. (D4) involves 
dpr{q,uj)/dq, some changes in the expression on the right-hand side of Eq. ( pl)l[ ) are necessary. To appreciate this, 
note that — rvQq)/dq — —rvQS{Lu — rv^q), from which, on account of the term {uj — rvQq)'^°~^ , one immediately 

observes that independent of the value of g, dpr{q,uj)/dq = oo. This would imply a peculiar situation where pr{q,uj), 
which is continuous almost everywhere, would be nowhere on the g-axis diffcrcntiable. As we shall see, this is wholly 
attributable to the fact that the expressions on the right-hand side of Eq. ( plD which are multiplied by distributions 
Q{uj =F rvoq) and Q{—uj ± rvoq) do not qualify as test [unctions (Gelfand and Shilov 1964). For simplicity, in the 
following we shall explicitly deal with the case corresponding to r = We write 

= 0(w + voq)Q{uj - voq)(p{q,ijj) + (uj ^ -^^,9 ^ -q), (D6) 
■.= Ai 01 + Voq) 4)2{uj - Voq), (D7) 

01 (x) 1=7(70,2^2;), (D9) 

/A \7o-i / A \ 
0.(x) :=(-.) exp (--.). (DIO) 

Since < 70 < 1/2, it is evident that the problem at hand (i.e. dp+{q, uj)/dq = 00) has its origin in the singularity of 
4>2{x) at a; 0. Through apphcation of the prescription as specified in Fig. 3 (see caption to this Figure), we define 
G+^i,{q, z) and consider G+{q, z) as being obtained through limj^jo G+^y{q, z). We have 

J-00 Z - LO Jvoq+u Z - Ul 

0(g,«og-Kz.)[l + (^^~^;og)H ^ r«^+'^ c^[q,voq + y) ^^^^^ 



It is from this expression that the derivative with respect to q has to be taken; we then have 
dG+{q,Lo) / dq:^\m\^lo dG+^^{q,uo) / dq. After some algebra, one obtains the following final result 

dG+{q,z) _ H~q,voq) , , d<j>{~q, ~u^')ldq 

Vo ; r I cl^ ; 



dq z -f Voq 

^ , {dd^^ju' + voq)lduo'} 02(^' - ^Qg) 
+ AlVo / do; 

J voq ^ - ^ 

d ( (I)i{lo' + Voq) ' 



+ Mvo ^^'h[^'-Voq)^{^-^^^^^V (D12) 

J voq OUJ (. Z — UJ J 

In arriving at this expression, we have made use of the properties d(j)i,2{(^' ± voq)/dq = ±vo d4>i,2{'^' ± voq)/duj' . It 
is evident that the integrands of the last two integrals on the right-hand side of Eq. ( D12| ) are singular (owing to 



02 ('^' — Voq)) at the lower limits of the integration boundaries; however since < 70 < 1/2, these singularities are 
integrable and therefore dG+{q, z)/dq as presented in Eq. ( D12 ) is a well-defined function of q and z. 



We are now in a position to draw the important conclusion that since (j){—q,voq) = ■Ai0i(O)02(2fog) oc <l"'"~^, 
dG+{q, z)/dq diverges as g | 0, for all z. This is in full conformity with our statement in Appendix C that, when 
dG+{q, z)/dq is unbounded at a particular value of g, it is unbounded at that g for all z. 



2. Self-energy and its asymptotic forms close to the Fermi points and the Fermi energy 



Now we proceed with evaluating the asymptotic behaviour of I]r(g,w) for w — > 0. We consider two cases, corre- 
sponding to g = and g 7^ 0. We demonstrate that, as 70 approaches zero, the self-energy of the Luttinger (1963) 
model becomes vanishingly small in a manner that is specific to non-interacting Fermi systems. 

Since we are interested in cases where < 70 < 1/2, for the determination of the leading asymptotic contribution 
to Gr{q,Lo), and thus to Sj.(g,aj), for g ^ (i.e. k kp) and uj ^ Q (i.e. e — > ep) we need to merely take account of 
the w'-integration in Eq. (D4) restricted to the interval [—A, A], where A is some finite constant satisfying |a;| < A 
(see Appendix C). Thus we define (cf. Eq. (D4)) 
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-A Z-U' 



Below we shall deal with ci^'' (q, lo) 



(D13) 



a. The q — case 



From the expression in Eq. ( |Dl| ) it can readily be deduced that (see Eqs. ([D6|)- (|D10[) ) 

Pr{q = 0,w) - Aya+i \uj\^'^°^\ for Lo^O. 
Making use of the general result < A) 

TT 



(DM) 



V [ duj' 
Jo 



= — tan — \u) 



~uj^ 2 



2 / ' ' a -I 



2F1 \ 1 



TT /7rcr\ „ 1 A'^-l A'^-3 
-tan(^— 1 u"^^ 



2 y 



(7—1 (7 — 3 
CT > -1, (T 7^ 1,3, 



(D15) 



where 2^1(0, 6; c;z) stands for the Gauss Hypergeometric function (Abramowitz and Stegun 1972, p. 556), for < 
7o < 1 and lo-"! < A we obtain 



G(^)((? = 0,..) ^ -7r^,„+i{ tan (!L(^2|_J)) + ,} sgn(c.)|^p^»-i 



(D16) 



Evidently, for < 70 < 1/2, the right-hand side of this expression diverges as w — > 0. As a consequence, for 
< 7o < 1/2, the possibly non-vanishing constant {Gr{q = 0,a;) — Gi^\q — 0,cj)} is asymptotically irrelevant for 
cij — > and we can write 



G,,((/ = 0,w)~^^^„+i{cot(^7o)-i}sgn(c^)|c^|2T«-i, lo ^ Q (0 < 70 < 1/2). 
From the Dyson equation we have 



(D17) 
(D18) 



S,. (g, w) = Gq.^ {q,uj) - G^'^{q,uj). 
Since Go;r(97^) = l/[^ ^ ''^o'rq + isgn(wor(j)] (see footnote |65|), it follows that 

G^^.{q,Lo) ^ LO - varq, (D19) 

and thus Gq!^((/ = 0,cij) = lo. For < 70 < 1, the leading-order contribution to Yir{q = 0,cij) is therefore entirely due 
to G:^^{q = b, w) on the right-hand side of Eq. (|D18|) . From Eq. ( |D16[ ) we thus obtain (0 < 70 < 1/2) 

1 



G;1((Z = 0,u;) 



-^rl cot(7r7o) + i| sgn(tj)|w|^" 



270 



7r^-y(,+i(cot^(7r7o) -I- l) 

As a consequence of our above arguments, for < 70 < 1/2 we have (for 70 in this range, cot(7r7o) > 0) 

s;((,.o,c.)^-^^^^i(^^^ 

7r^^o+i(cot^(7r7o) -|- 1) 
7ryl^„+i(cot^(7r7o) -f 11 



1 1-270 



(D20) 

(D21) 
(D22) 



We note that Sr('? = 0,0-1 = 0) = (see Eq. ( D17D and text preceding it), which in view of our considerations in 
Appendix B amounts to the fact that in the Luttinger model Tj^^{q = 0) is entirely cancelled by correlation effects. 
It is further interesting to note that 

S;'((? = 0,u;) 



tan(7r7o), a; — + 0, 



(D23) 



E;(g = 0,0;) - E^(g = 0,o; = 0) 
which is seen to approach zero (like ^ 7r7o) for 70 i 0. As we have mentioned in § IV. C, 70 = corresponds to the 



non-interacting system, so that the expressions presented in Eqs. (D21) and (D22) reproduce the non-interacting limit 
in a continuous way. 
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6. The q ^ case 



Here we confine our considerations to the case corresponding to r = +. From Eq. ( Dll ) (taking the limit J, 0) it 
readily follows that for < 70 < 1 the leading asymptotic contribution to G+{q,uj), as w 0, is due to 



(D24) 



voq 



For z = voq + lj + irj, rj I 0, with < w < A, making use of 7(70, A(a;' + 2voq)/[2vo]) — 7(70, Ag) + 0(10'), 
exp(— Acj'/[2wo]) = 1 + 0{uj'), while employing the standard result {0 < ui < A) 



r 



= — TT COt(7r(7)ti>'^ 



cr- 1 



2Fi(l,l-a;2 



1 - cr 2 - (7 
(7 > 0, cr ^ 1,2, 



60 + 0(^2-)^ 



we obtain (0 < w < A, < 70 < 1 and g ^ 0) 



,(A) 



{q,uj) - G'+(g,w) - 7r7(7o,Ag) | cot(7r7o) -iji 



,70-1 



W i 0. 



(D25) 



(D26) 



Note that the exponent 70 — 1, as compared with 270 — 1 for the case corresponding q — Q, implies that even for 70 
in the open interval (0, 1) (as opposed to (0, 1/2)), G+{q 7^ 0, w) is to leading order asymptotically identical with the 
corresponding auxiliary function Q^'^^q ^Q,lo) as uj [ 0; recall that the leading asymptotic terms of Gr^\q = 0,a;) 
and Gr{q — 0,a;) for w ^ (see Eqs. (D16) and (D17)) are identical only when < 70 < 1/2. Further, similar 
to the case corresponding to g = 0, the divergence of G'+(g, w) for < 70 < 1, as j 0, implies that the possibly 
non-vanishing constant {G'_|_((7 7^ 0,a;) — Q')^^ {q 7^ 0,a;)} is asymptotically irrelevant. 
Proceeding as in the case of 9 = 0, for w J, we obtain (0 < 70 < 1) 



G-^\qM 



7^7 (70, Ag)^^o ( cot2(7r7o) + 



■jyl cot(7r7o) + i^u?' 



(D27) 



As a consequence of our above arguments, for < 70 < 1 we have (for 70 in this range, cot(7r7o) takes on both positive 
as well as negative values) 



S'+(<Z,c^) 



s^(g,^) 



cot(7r7o) 



777(70, Ag)^^o(cot2(7r7o) + l) 
-1 

777(70, Ag)^^o(cot^ (7770) + 1) 



,1-70 



,1-70 



(D28) 
(D29) 



We point out that I]+((j,a; = 0) = (see Eq. ( |D17| ) and text preceding it), which in view of our considerations in 
Appendix B leads to the conclusion that, in the Luttinger model Y?^^ (jq) is entirely cancelled by correlation effects. 
As in the case corresponding to g = 0, here we have (compare with Eq. ( |D23| )) 



S^(<Z,^) 



E;(g,co)-E+(q,L^ = 0) 



tan(7r7o), i 0, 



(D30) 



which approaches zero (like ^ 7770) for 70 i 0; this result is consistent with the fact that 70 = corresponds to the 
non-interacting system. 

We further point out that since d^(a,z)ldz ~ exp(— z)z'^^^, from Eqs. (D28) and ( D29| ) we obtain that 
dYij^-{q,bj) / dq oc which diverges as g — > for 70 < 1. This result is in full conformity with our finding 

above (see the last paragraph in § D.l, following Eq. (DIS)). 

In closing this Appendix, we point out that the apparent difference in the asymptotic behaviours of T,.r{q,Lo) 
(specifically concerning r = + explicitly considered here) for q ^ and q = 0, is not specific to the Luttinger 
model. It is well known (Hodges, Smith and Wilkins 1971, Bloom 1975, Fujimoto 1990, Fukuyama, Narikiyo and 
Hasegawa 1991, Fukuyama, Hasegawa and Narikiyo 1991) that for isotropic systems of interacting electrons in two 
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spatial dimensions, under the assumption that they are Fermi liquids, while £)~'(e — for k ^ kp, we have 

S"(fcF;e)~'(£-£F)^ln|£-£F|. □ 
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k 




FIG. 1. A second-order skeleton self-energy diagram in terms of the bare particle-particle interaction v (broken lines) em- 
ployed by Luttinger (1961). The solid lines stand for either Go or G, the single-particle Green function corresponding to the 

non-interacting and interacting systems, respectively. The external and internal wave- vectors, k and ki, k2 and ka, respectively, 
are shown. Because of conservation of momentum, one of the internal wave- vectors, say ka, can be eliminated; ka = ki -|-k2 — k. 




FIG. 2. Contour C of integration employed in the proof of the Migdal theorem concerning the discontinuity of the momentum 
distribution function n{k) aX k = kp and its relation to the quasi-particle weight on the Fermi surface, Z^^. The shaded parts of 
the real axis indicate branch cuts of G(A;; z); the unshaded paxt of this axis centred around the 'chemical potential' n represents 
{fj,N,fiN+i) with fj,N and fiN+i two branch points of G{k;z) as well as T,{k;z). 
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FIG. 3. The distribution (as opposed to function) f{ij):—0(ij — u}o)g(ij) and a specific function fviuj), corresponding to a 
finite positive value of i/, in terms of which we choose fiuj) to be defined (this choice is not unique); /(ci;)~ lim^jo fv{ijj)- We 
have the following /i/(a;) = 0, for lo < — v; fu{ijj) = g{LOo + + (w — wo)/^; ioi uJo — v < co < u>o', fviijj) = g{>jJo + i'), for 
u)o < ij) < u)o + v and /i/(w) = gi^jj), for ui > u)o + v. The form of f{u)) is reminiscent of p+(Q,a;) pertaining to the Luttinger 
model for which wo coincides with voq and p+(g, w) oc (a; — voqy°~^ for w J. voq (here < 70 < 1/2). 
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